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Galaxy evolution in Hickson compact groups: The role of ram 
pressure stripping and strangulation 
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ABSTRACT 

Galaxies in compact groups tend to be deficient in neutral hydrogen compared to isolated 
galaxies of similar optical properties. In order to investigate the role played by a hot intragroup 
medium (IGM) for the removal and destruction of Hi in these systems, we have performed a 
Chandra and XMM-Newton study of eight of the most Hi deficient Hickson compact groups. 
Diffuse X-ray emission associated with an IGM is detected in four of the groups, suggesting 
that galaxy-IGM interactions are not the dominant mechanism driving cold gas out of the 
group members. No clear evidence is seen for any of the members being currently stripped of 
any hot gas, nor for galaxies to show enhanced nuclear X-ray activity in the X-ray bright or 
most Hi deficient groups. Combining the inferred IGM distributions with analytical models 
of representative disc galaxies orbiting within each group, we estimate the Hi mass loss due 
to ram pressure and viscous stripping. While these processes are generally insufficient to 
explain observed Hi deficiencies, they could still be important for Hi removal in the X-ray 
bright groups, potentially removing more than half of the ISM in the X-ray bright HCG97. 
Ram pressure may also have facilitated strangulation through the removal of galactic coronal 
gas. In X-ray undetected groups, tidal interactions could be playing a prominent role, but it 
remains an open question whether they can fully account for the observed Hi deficiencies. 

Key words: galaxies: evolution — galaxies: interactions — galaxies: ISM — X-rays: galaxies 
— X-rays: galaxies: clusters 



1 INTRODUCTION 



The origin of the galaxy morphology-density relation is still one 
of the most important unsolved problems in astrophysics. Not only 
are spiral galaxies less common within dense cluster environments, 
but those which are present tend to be deficient in Hi, and this de- 
ficiency itself correlates wi th projected local galaxy density (e.g., 
iGiovanelli & Havneslll985h . The mechanisms responsible for the 
changes in the morphology and gas content of galaxies are unclear, 
with gas stripping, tidal shocks, and galaxy interactions and merg- 
ers all contenders. 

Traditionally, clusters of galaxies have represented the envi- 
ronment of choice for attempts to unravel the nature of the rele- 
vant processes. Recent results, however, strongly suggest that the 
origin of the environmental modification of galaxies which under- 
pin the morphology-density relation lies, not in the cores of galaxy 
clusters, but in smaller groups and cluster outskirts. For exam- 
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pie, spectroscopic studies of the effe cts of the cluster environment 
on galaxies (e.g., iLewis et alj|2002h show that the suppression of 
star formation takes place in cluster outskirts rather than in the 
core, and is modulated by local galaxy density. Moreover, X-ray 
bright groups sh ow a morphology-density relation stronger than 
that of clusters (Hels don & Ponmarj|2003h . a result adding to the 
accumulating evidence that cluster galaxies have often been 'pre- 
processe d' in groups prior t o their assembly into larger systems 
(see, e.g JCortese et alj2006t) . In addition, it is becoming clear that 
processes once thought to be ex clusive to the cluster environment, 
such as ram pressure stri pping dRasmussen. Ponman & Mulchaey 
2006) and strangulation JKawata & Mulch aey 2008), may play a 
role also in much smaller systems. In order to elucidate the ori- 
gin of the morphology-density relation, it is therefore necessary to 
study the processes acting on galaxies within groups. 



The compact groups in the catalogue of Hicksor] dl982l) of- 
fer particularly interesting opportunities in this regard. Many of 
these groups are spiral-rich, but their galaxy population is, on av- 
erage, deficient in Hi by a factor of ~ 2 compared to loose groups 
(Williams & Roocj |l987l) . Furthermore, some of these groups have 
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exceptionally compact galaxy configurations, and so may repre- 
sent pre-virialisation systems close to maximum collapse, in which 
galaxies are suffering strong environmental modification but have 
yet to be converted into early-types. While recent work has iden- 
tified tidal interactions and mergers as playing an important role 
in the morphological transformation o f spirals in some compact 
groups jCoziol & Plauchu-Fravrj|2007h . it is still unclear what is 
causing the observed Hi deficiencies and to what extent this is re- 
lated to the processes modifying the stellar component of the group 

me mbers. 

IVerdes -Montene gro et alj d200lh presented a detailed study of 
the Hi content of a sample of 72 Hickson compact groups (HCGs), 
with integrated Hi masses from single dish measurements, and de- 
tailed Very Large Array (VLA) mapping of a subsample. Defining 
Hi deficiency Am as 

Am = logMm.prcd - logMm.obs, (1) 

where Mm.obs is me observed Hi mass of the group galaxies and 
Mm.pred is that predicted f or isolated galaxies of simil ar morphol- 
ogy and optical luminosity (Haynes & Giovanelli 1984), this work 
confirmed the deficiency in HCGs, and allowed a search for corre- 
lations between deficiency and other group properties. One of the 
strongest relationships found was that with detectable intergalactic 
X-ray emission; almost half of the groups with significant Hi de- 
ficiency showed diffuse intragroup X-ray emission in the ROSAT 
survey of Ponm anet alj i 19961 hereafter P96). More recently, sim- 
ilar result s have been reported also for groups outside Hickson's 
catalogue ( Sengupta & Balasubramanvam 2006). 

The increased prevalence of X-ray detected systems among 
Hi deficient (compact) groups may suggest a picture whereby 
Hi is stripped from spiral galaxies within virialising groups, and 
then destroyed due to heating by a surrounding hot intragroup 
medium (IGM). However , for m ost of the compact groups in the 
Verdes-Montenegro et al. I fcooih sample, ROSAT data were either 
not available or of insufficient quality to permit any detailed study 
of the relationship between the hot and cold gas. With only very 
shallow ROSAT All-Sky Survey data at their dispo sal for a number 
of these svstems. lVerdes-Montenegroe t al. (2001) could not estab- 
lish the amount and properties of any hot IGM in these groups. 

In order to explore the processes destroying Hi in these sys- 
tems, we have therefore embarked on a programme to obtain high 
quality X-ray and radio data for the most Hi deficient compact 
groups, adding in existing archival X-ray data wherever relevant. 
Since ram pressure stripping can only operate where a significant 
IGM is present, while tidal stripping of gas requires only that galax- 
ies be interacting, a key discriminator for the mechanism of Hi re- 
moval from galaxies is whether or not there is a correlation between 
Hi deficiency and the properties of a hot IGM, and in particular 
whether a hot IGM is present in the highly Hi deficient systems. 
Clarifying these issues represents the goal of the present study, in 
an attempt to shed light on the role played by galaxy-IGM inter- 
actions in destroying Hi and establishing the morphology-density 
relation in these compact groups. The focus of this paper thus 
rests mainly on the X-ray properties of any hot intergalactic gas 
within the groups, while a forthcoming companion paper (Verdes- 
Montenegro et al., in prep.) will discuss in more detail the Hi and 
radio continuum emission in the groups, the detailed relationship 
between X-ray and Hi morphology, and the multi-wavelength prop- 
erties of the individual group galaxies. 

In Sections [2] and [3] we outline the sample selection and X- 
ray data analysis, respectively. Section [4] presents the results for 
the X-ray properties of the hot gas and galaxies within each group, 



and compares the derived IGM properties to the observed Hi defi- 
ciencies. In Section[5]we construct an analytical model of a repre- 
sentative late-type galaxy orbiting within the derived gravitational 
potential of each X-ray bright group, allowing us to evaluate the 
importance of ram pressure and viscous stripping for Hi removal. 
The results are discussed in Section[6]and summarized along with 
the main conclusions in Section|7] 

A Hubble constant of 73 km s~ Mpc" 1 is assumed through- 
out. Unless otherwise stated, all errors are quoted at the 68 per cent 
confidence level. 



2 GROUP SAMPLE AND OBSERVATIONS 

Our broad aim is to establish the relationship between the hot and 
cold gas in Hi deficient groups, and explore the processes of gas 
removal operati ng within them. Our sample is t herefore based on 
that studied by IVerdes-Mo ntenegro et all J200 lh . from which we 
selected all HCGs which are highly deficient in Hi (Am > 0.5) 
based on Very Large Array (VLA) Hi data, contain four or more 
group galaxies, and lie at a distance D < 100 Mpc. This yielded an 
initial sample of 1 1 groups, which will be discussed in its entirety 
in Verdes-Montenegro et al. (in prep.). VLA Hi mapping has been 
completed for all 11 systems, along with follow-up Green Bank 
Telescope (GBT) observations for the groups discussed in this pa- 
per. Compared to the VLA, the single-dish GBT is more sensitive 
to extended, smoothly distributed emission that may otherwise be 
filtered out by the VLA interferometer. Including the GBT data thus 
provides a more complete census of the Hi content in the groups. 
For HCG 48, included in the initial sample, the additional gas de- 
tected by GBT indicates that this group is not Hi deficient after 
all. Consequently, this group was omitted from the sample for the 
purpose of the present study. 

Of the remaining systems, the eight included in this paper are 
those for which Chandra or XMM-Newton X-ray data are currently 
available. Four of these eight groups a re exceptionall y compact, 
with major galaxies (as catalogued by lfficksoni ri982) contained 
within a circle of radius ~ 2 arcmin, and thus requiring Chan- 
dra data to resolve their X-ray structure in the crucial region. We 
note that these groups are not generally mature, X-ray bright sys- 
tems dominated by early-type galaxies. Such groups may contain 
little Hi, but the Hi content of their galaxies as predicted from the 
lHavnes & Giovanellil ( fl984l) results mentioned above is also small 
(zero for ellipticals and modest for lenticulars). Such groups are 
therefore not usually Hi deficient according to our definition in 
equation QJ. Our chosen systems have Hi content far below that 
expected for their galaxy contents, and are therefore those in which 
the processes which destroy Hi should be in active, or very recent, 
operation. 

The observation log for the X-ray data considered in this paper 
is presented in Table Q] detailing the pointing coordinates for each 
observation (for archival observations not necessarily identical to 
the optical group centre), group distance D, the observing instru- 
ment, date, and mode, along with cleaned exposure times t oxp for 
each camera, and the Gala ctic absorbing column density Nh from 
iDickev & Lockmanl J 1990) as adopted in the X-ray spectral analy- 
sis. 
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Table 1. Log of available X-ray observations. Group luminosity distances D for the adopted value of Hq were taken from 
the NASA/IPAC Extragalactic Database (NED). Column 7 specifies the frame mode (full frame/extended full frame) and 
optical blocking filter for XMM, and ACIS CCD aimpoint and telemetry mode (Faint/Very Faint) for Chandra. 



Group 


RA 


Dec 


D 


Chondral 


Obs. date 


Obs. mode 


^cxp 






(J2000) 


(J2000) 


(Mpc) 


XMM 


(yyyy-mm-dd) 




(ks) 


(10 20 cm" 2 ) 


HCG 7 


no i ^ s 






VMM nn 
ALVllVl pil 


9004 1 ? ")(\ 


FP Thin 

rr i n in 




2 24 










VMM m 1 

JllVllVI 1111 




FF Thin 


35 1 












XMM ml 




FFThin 


35.9 




HCG15 


02 07 39.0 


+02 08 18.0 


92 


XMM pn 


2002-01-10 


EFF Thin 


23.3 


3.20 










XMM ml 




FFThin 


31.0 












XMM m2 




FFThin 


31.1 




HCG 30 


04 36 28.3 


-02 50 02.9 


63 


Chandra 


2006-02-07 


ACIS-S VF 


29.0 


5.08 


HCG 37 


09 13 35.3 


+30 01 25.5 


97 


Chandra 


2005-01-13 


ACIS-S VF 


17.5 


2.30 


HCG 40 


09 38 56.7 


-04 50 55.9 


98 


Chandra 


2005-01-29 


ACIS-S VF 


31.8 


3.64 












2005-01-29 


ACIS-S VF 


14.6 




HCG 44 


10 18 02.5 


+21 48 50.7 


23 


Chandra 


2002-03-14 


ACIS-S F 


19.4 


2.16 




10 17 38.0 


+2141 17.0 




XMM pn 


2001-05-07 


EFF Thick 


8.1 












XMM ml 




FFThin 


13.4 












XMM ml 




FFThin 


13.1 




HCG 97 


23 47 25.4 


-02 19 45.5 


86 


Chandra 


2005-01-14 


ACIS-S VF 


36.2 


3.65 


HCG 100 


00 01 25.9 


+13 06 30.4 


69 


Chandra 


2006-06-12 


ACIS-I VF 


26.1 


4.40 












2007-01-24 


ACIS-I VF 


16.1 





3 DATA REDUCTION AND ANALYSIS 
3.1 XMM-Newton data 

The XMM data were analysed using XMMSAS v6.5.0, and cali- 
brated event lists were generated with the 'emchain' and 'epchain' 
tasks. Event files were filtered using standard quality flags, while 
retaining only patterns ^ 4 for pn and 12 for MOS. Screening 
for background flares was first performed in the 10-15 keV band 
for MOS and 12-14 keV for pn. Following an initial removal of 
obvious large flares, a 3a clipping of the resulting lightcurve was 
applied. Point sources were then identified by combining the re- 
sults of a sliding-cell search ('eboxdetect') and a maximum likeli- 
hood point spread function fitting ('emldetect'), both performed in 
five separate energy bands spanning the total range 0.3-12 keV. In 
order to filter out any remaining soft protons in the data, a second 
lightcurve (3<r) cleaning was then done in the 0.4—10 keV band, 
within a 9-12 arcmin annulus which excluded the detected point 
sources. Closed-filter data from the calibratio n database and blank- 
sky background data ( Read & Ponman 200^) for the appropriate 
observing mode were filtered similarly to source data, and screened 
so as to contain only periods with count rates within lcr from the 
mean of the source data. All point sources were excised out to at 
least 25 arcsec in spectral analysis. 

To aid the search for diffuse X-ray emission within the groups, 
smoothed exposure-corrected images were produced, with back- 
ground maps generated from blank-sky data. We allowed for a 
differing contribution from the non-vignetted particle background 
component in source- and blank-sky data by adopting the following 
approach. First an EPIC mosaic image was smoothed adaptively 
(3<7-5<t significance range), and the particle background was sub- 
tracted. The latter was estimated from closed-filter data which were 
scaled to match source data count rates in the image energy band 
in regions outside the field of view, and then smoothed at the same 
spatial scales as the source data. The resulting photon image in- 
cludes the X-ray background at the source position. To remove this 
component, a particle-subtracted blank-sky image was produced in 
a similar way, and scaled to match 0.3-2 keV source count rates 
in a point-source-excised 10-12 arcmin annulus assumed to be 



free of IGM emission (this assumption is clearly justified in all 
cases, as will be shown). This image was then subtracted from the 
corresponding source image, and the result was finally exposure- 
corrected. 



3.2 Chandra data 

For all Chandra data sets, calibrated event lists were regenerated 
using CIAO v3.3. For Very Faint mode observations, the standard 
additional background screening was carried out. Bad pixels were 
screened out using the bad pixel map provided by the pipeline, 
and remaining events were grade filtered, excluding ASCA grades 
1, 5, and 7. Periods of high background were filtered using 3<r 
clipping of full-chip lightcurves, binned in time bins of length 
259. 8-s and extracted in off-source regions in the 2.5-7 keV band 
for back-illuminated chips and 0.3-12 keV for front-illuminated 
chips. Blank-sky background data from the calibration database 
were screened and filtered as for source data, and reprojected to 
match the aspect solution of the latter. Point source searches were 
carried out with the CIAO task 'wavdetect' using a range of scales 
and detection thresholds, and results were combined. Source ex- 
tents were quantified using the 4cr detection ellipses from 'wavde- 
tect', and these regions were masked out in all spectral analysis. 

In order to produce smoothed images as for the XMM data, 
background maps were generated using blank-sky data, and scaled 
to match source count rates for each CCD. This scaling employed 
either the full-chip 10-12 keV count rates, with point sources ex- 
cluded, or, where possible, count rates in the image energy band 
within source-free regions on the relevant CCD. The background 
maps were then smoothed to the same spatial scales as the source 
data and subtracted from the latter. The resulting images were 
finally exposure corrected using similarly smoothed, spectrally 
weighted exposure maps, with weights derived from spectral fits to 
the integrated diffuse emission (where possible, otherwise the ex- 
posure maps were weighted by a T — 1 keV, Z = 0.3 Zq thermal 
plasma model). 
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3.3 Spatial and spectral analysis 

While the smoothed X-ray images described above are useful in 
terms of establishing the presence and rough morphology of any 
intragroup medium, we emphasize that they were used for illus- 
trative purposes only and not for quantitative analyses. Where 
IGM emission was not immediately obvious from these images, 
we performed an additional source detection procedure based on 
Voronoi tessellation and percolation ('vtpdetect' in CIAO), which 
can be useful for detecting extended, low-surface brightness emis- 
sion missed by our standard detection algorithms. To reduce the 
fraction of spurious detections, a minimum of 50 net counts were 
required for a source to be considered real. 

As a second step in the search for group-scale diffuse emis- 
sion, we also extracted exposure-corrected 0.3-2 keV surface 
brightness profiles from the unsmoothed data, with all detected 
point-like and extended galactic sources masked out. The profiles 
were extracted from t he optica l group centre defined by the prin- 
cipal members in the iHicksonl dl982h catalogue. For XMM data, 
we used an EPIC mosaic image for this purpose, with the particle 
background removed using the method described above. The es- 
timated particle level shows a typical standard error of the mean 
of « 5 per cent, which should be representative of the uncertainty 
associated with particle subtraction if the ratio of particle events 
inside and outside the field of view is similar to that in the closed- 
filter data. We used the blank-sky background data to confirm this 
assumption (since these have very little source contamination), but 
have conservatively added a 10 per cent error in quadrature to our 
XMM surface brightness errors, to allow for any residual systematic 
uncertainties associated with the particle subtraction. 

For the spectral analysis of any extended emission, X-ray 
spectra were accumulated in energy bins of at least 20 net counts, 
and fitted in XSPEC vll.3 assuming an APEC thermal plasm a 
model with the solar abundance table of lGrevesse & Sauvall Jl998t> . 
XMM background spectra were extracted by means o f the common 
'double-subtraction' technique dArnaud et ai1l2002h . using blank- 
sky background data for the on-chip background, and a large-radius 
(10-12 arcmin) annulus for determining the local soft X-ray back- 
ground. Owing to the smaller field of view, a similar approach was 
not generally possible or desirable for the Chandra observations 
where source emission may completely fill the CCD under inves- 
tigation. The extraction of Chandra background data products are 
therefore described individually for each group in the next Section. 

Surface brightness profiles of the X-ray detected groups were 
extracted from the peak of the diffuse X-ray emission when clearly 
identifiable (in HCG 37 and 97) and from the centroid otherwise 
(HCG 15 and 40). The profiles were fitted with standard /3-models 
for conversion into IGM density profiles under the assumption of 
isothermality. Since X-ray emissivity is very nearly independent of 
temperature for a T ~ 0.5-1 keV pla sma of the relevant metal - 
licities (see [Sutherland & Dopitalll993l) , this approach is entirely 
adequate for our purposes, where the uncertainties of our analysis 
will ultimately be dominated by those related to the modelling of 
the impact of the hot gas on the galaxies. The density profiles were 
normalized using the spectral normalization A from XSPEC, 

10~ 14 f 

A= 4nD>{l + z)' y^Hdl/cm-5, (2) 

where D is the assumed group distance, and n e and nu are the 
number densities of electrons and hydrogen atoms, respectively. 
The integral represents the fitted emission integral I c over the cov- 
ered volume V, assumed to be spherically symmetric. Total gas 



masses within the region of interest were derived by simple vol- 
ume integration of n e (r)/i e m p , where m p is the proton mass and 
we have assumed n e /riH = 1.165 and /j, e — 1.15, appropriate for 
a fully ionized Z = 0.3 plasma at the relevant temperatures 
dSufherland & Dopitall 19931) . 

Since we have no knowledge of the density distribution of 
any hot gas in the X-ray undetected groups, we have only derived 
constraints on their mean electron density (n e ) within the region 
considered, effectively assuming a uniform IGM distribution. The 
advantage of this approach is that it provides very conservative up- 
per limits to the IGM masses (and mean ram pressures) within the 
relevant region. The derived limits to (n e ) and Migm for these 
groups were obtained from the IGM count rate limits and thus de- 
pend on the depth of the X-ray data. These count rate limits were 
established on the basis of the exposure-corrected (and, in the case 
of XMM, particle-subtracted) images, by comparing the emission 
level within a region centred on the optical group centre with that in 
a surrounding annulus. The physical extent of the region of interest 
was thus constrained by the need to evaluate the background locally 
from our data, and varies from 50-150 kpc among the groups, as 
detailed in the discussion of individual groups below. The derived 
constraints on IGM count rates were translat ed into constraints on 
Z e , ass uming the cooling function A(T, Z) of lSutherland & Dopital 
d 19931) and an IGM t emper ature taken from the cry-Tx relation of 
lOsmond & Ponmart Jiooj) , 

log (j v = (1.15 ±0.26) log T x + 2.60 ±0.03, (3) 

with galaxy velocity dispersions ov in km s _1 and Tx in keV. Er- 
rors on T were derived from the dispersion of this re lation, with cry 
taken from P96 for HCG 7, 3 0, 37, 44, and 100, fromlMahdavi et all 
d2005l) for HCG 97, and from lOsmond & Ponmanl d2004l) for theie- 
mainder. The resulting temperature range was then used to estimate 
3<r upper limits on (n e ) ~ (Ic/V) 1 ^ 2 inside the assumed spheri- 
cal volume V for any subsolar metallicity Z. The assumption of a 
uniform IGM in these groups implies that we can identify the up- 
per limit to the central IGM density no with (n e ) and constrain the 
IGM masses by simply multiplying (n e ) and V . 

In order to briefly investigate the level of nuclear X-ray activ- 
ity among individual group members, 0.3-2 keV count rates of all 
galactic central point sources were also extracted, adopting extrac- 
tion regions of 2 and 15 arcsec radius for Chandra and XMM data, 
respectively. In the majority of cases, photon statistics were insuf- 
ficient to allow robust spectral fitting for individual sources. For 
consistency, all point source count rates were therefore converted 
to luminosities assuming a power-law spectrum of photon index 
T = 1.7, absorbed by the Galactic value of A?h- The associated 
uncertainties were derived from the Poisson errors on the photon 
count rates. 



4 RESULTS 

In this section we discuss the results obtained for the IGM in each 
group and for the X-ray emission associated with individual group 
galaxies. Figure Q] shows contours of the smoothed, background- 
subtracted X-ray emission of each group overlayed on Digitized 
Sky Survey (DSS) images. Diffuse X-ray emission associated with 
an intragroup medium is detected in four of the eight groups, as de- 
scribed for each group individually below. For the remaining four, 
we do not detect any extended group emission, neither inside a 
given physical radius from the optical group centre when compared 
to the emission level in surrounding regions, nor on the basis of the 
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Figure 1. Con tours of a daptively smoothed 0.3-2 keV emission overlayed on DSS images for all groups, with the principal group members labelled following 
the notation of lHicksonl Jl982l) . Where relevant, dashed squares outline the coverage of the Chandra/ ACIS CCD's. For HCG44, dark (black) contours outline 
the Chandra emission, and lighter (green) those of an overlapping XMM pointing. 



Voronoi source detection procedure. To corroborate these results, 
Figures[2]and[3]show the derived surface brightness profiles for the 
groups with and without detectable intragroup emission, respec- 
tively. 

Table [2] summarizes the observed Hi and X-ray properties of 
the groups, along with the adopted velocity dispersions from optical 
spectroscopy. Hi deficiencies in the Table are from our GBT mea- 
surements (Borthakur et al., in prep.), exc ept for HCG 44, for which 
we h ave adopted the older VLA value iVerdes-Montenegro et al.l 
l200ll) due to the GBT beam size only covering the central region 
of this relatively nearby system. The listed Hi deficiencies are based 
on the integrated Hi mass within the circular region covered by the 
radio data (rm in the Table). In many groups, a significant fraction 
of the detected Hi is located outside the optical extent of individ- 
ual galaxies (i.e. is intergalactic) and cannot be clearly assigned to 
any individual group member. The derived values of Am should 



therefore generally be viewed as an average for the galaxies within 
the GBT or VLA field. As the fractional 1-cr uncertainty on the 
measured Hi masses from our GBT data is less than 1 per cent for 
all groups, uncertainties on the listed deficiencies are dominated by 
those related to the predicted logarithmic Hi mass, which we have 
taken to be 0.2, adopt i ng the standard estimate of error provided by 
lHavnes & Giovanellil < fl984T) . 

For reference, 0.3-2 keV X-ray luminosities are also listed in 
Table [2] corrected for Galactic absorption, and derived within the 
region employed for the spectral analysis unless otherwise speci- 
fied in the subsection for the relevant group. The listed central hot 
IGM densities (or upper limits to the mean densities for the X- 
ray undetected groups) were computed as outlined in Section [331 
IGM masses in the Table were derived within the same region as 
the Hi deficiencies (i.e. within rui given in the Table), to enable a 
direct comparison between the two. Note, as indicated above, that 
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Table 2. Summary of derived group properties. Except where indicated, Hi deficiencies Ajji are from GBT 
data (Borthakur et al., in prep.), obtained inside a radius rjji. Column 7 lists the derived constraints on central 
hot IGM density, with upper limits for the X-ray undetected groups given at 3<r significance, and Column 8 
the corresponding hot IGM mass within rjji ■ 



Group 


A m 


(kpc) 


(tar 1 ) 


(keV) 


L X 

(10 40 erg s- 1 ) 


n 

(10~ 3 cm" 3 ) 


A^IGM 
(10 10 M ) 


HCG7 


0.60 


68 


95" 


0.3 ±0.1* 


< 0.5 


< 0.06 


< 0.26 


HCG15 


0.46 


113 


404 6 


n oq+0.09 


32 ± 2 


4.4 ±0.5 


7.4 ±0.8 


HCG30 


1.37 


79 


72" 


0.2 ±0.1* 


< 1.3 


< 1.2 


< 8.3 


HCG37 


0.33 


119 


446" 


n sfi+ - 13 


16 ± 2 


38 ±6 


6.5 ± 1.0 


HCG40 


0.60 


121 


157 6 


59+ 11 


3.1 ±0.5 


1.1 ±0.2 


2.4 ±0.4 


HCG44 


0.69t 


101 


145" 


0.4 ±0.1* 


< 3.6 


< 0.1 


< 1.4 


HCG97 
HCG 100 


0.35 
0.27 


106 
86 


383 c 
100" 


Q7+ 014 

U - M '-0.12 

0.3 ±0.1* 
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Figure 2. 0.3-2 keV exposure-corrected surface brightness profiles of groups with detectable diffuse emission, along with best-fitting /3-models (solid lines). 
Horizontal dotted lines mark the estimated background level in each case. For each plot, the bottom panel shows fit residuals relative to the best-fitting model. 
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Figure 3. As Fig. [2] but for the groups without detectable diffuse emission. 



the upper limits to Migm for the X-ray undetected groups conser- 
vatively assume a uniform IGM distribution. If instead assuming a 
standard /3-model for the hot gas in these groups, with central den- 
sity equal to the inferred mean value (n e ) and with typical group 
values of, e.g., (3 = 0.5 and r c = 20 kpc, the derived IGM mass 
limits would be reduced by factors of 4-7. 



4.1 X-ray properties of the intragroup medium 

4.1.1 HCG 7 

This group remained X-ray undetected in shallow ROSAT All-Sky 
Survey (RASS) data. Despite the XMM data of this target represent- 
ing the deepest X-ray observation within our sample, no diffuse X- 
ray emission is detected in the group when comparing the emission 
level of the exposure-corrected and particle-subtracted 0.3-2 keV 
mosaic image inside r = 9 arcmin (r w 140 kpc) with that mea- 
sured in a surrounding annulus. This is corroborated by the derived 
surface brightness profile shown in Fig. [3] Although this profile 
does seem to hint at a weak signal inside r w 2 arcmin, the com- 
bined signal inside this region is significant at less than 1.8a, is 
not picked up by 'vtpdetect', and is not discernible in the smoothed 
image presented in Fig.[T] Thus, we conservatively treat it as a non- 
detection. 



2.6t°;I keV 

0-411q'q2 for an assumed abundance of 0.3 so- 



In fact, no extended emission unassociated with individual 
galaxies is identified by 'vtpdetect' within the 9 arcmin radius, with 
the exception of the source visible in Fig. Q] roughly ~ 5 arcmin 
south of the optical group centre. There are no optical/infra-red 
counterparts to this X-ray source listed in NED within a 2 arcmin 
diameter, despite the proximity of the group (D ~ 54 Mpc). The 
emission is detected out to 2.2 arcmin from the X-ray centroid at 
3a above the local background, confirming that it is clearly ex- 
tended. A thermal plasma model fit to the spectrum extracted from 
the pn data within r = 1.5 arcmin of the centroid provides an ac- 
ceptable fit, with a reduced xt = 0-92 for 24 degrees of freedom 
(d.o.f.). This yields a best-fitting temperature T 
and redshift z 

lar. A simple power-law model with Galactic absorption yields 
r = 1.84 ± 0.10 for the power-law index, but the fit is not pre- 
ferred to a thermal model (red. x 2 = 1-07 for 25 d.o.f, i.e. a 
change of A^ 2 = 4.8). The estimated temperature can be com- 
pared to that expected for the IGM in HCG 7 on the basis of its 
very low galaxy veloci ty dispersion, ay = 95 km s -1 , for which 
the av-Tx relation of lOsmond & Ponmanl d2004h would suggest 
only T = 0.3 ± 0.1 keV. Combined with the redshift estimate of 
z ~ 0.4, this strongly suggests that this emission is not associ- 
ated with HCG 7 itself. The X-ray centroid also coincides to within 
10 arcsec with an NVSS source with a 1.4 GHz flux of 19.4 mjy 
(corresponding to 1 x 10 25 W Hz -1 at z = 0.41), so the X-ray 
emission is conceivably associated with a z ~ 0.4 background 
cluster harbouring a central radio-loud galaxy. 

For T anywhere in the range 0.2-0.4 keV and assuming any 
subsolar metallicity, our failure to detect IGM emission inside r « 
150 kpc translates into a 3(7 upper limit to the unabsorbed 0.3- 
2 keV luminosity inside this region of Lx < 5 x 10 39 erg s _1 , 
with a corresponding limit to the mean gas density of (n e ) < 6 x 
1 0~ 5 an" 3 . We not e that HCG 7 is included in the group catalogue 
of I Yang etai] ( l2007h . with a total group mass, as estimated from its 
optical properties, ranging from 2.4-5.0 xlO 12 Mq depending on 
the method assumed. This places HCG 7 at the very low-mass end 
of the group mass function, with a mass similar to that of the Local 
Group. Thus, it is perhaps not surprising that we fail to detect any 
IGM emission in this system. 



4.7.2 HCG 15 

This group has a relatively high velocity dispersion of ctv ~ 
400 km s _1 , and hot intragroup gas was already detected in pointed 
ROSAT observations (P96). Significant IGM emission is seen in the 
XMM data presented in Fig.Q] revealing a somewhat disturbed X- 
ray morphology. Despite this irregularity of the emission, an op- 
tically bright early-type galaxy is present roughly at centre of the 
X-ray emission, as is typical for fairly undisturbed X-ray bright 
groups. The facts that this galaxy is a lenticular rather than an ellip- 
tical, and that the IGM emission is not (yet) strongly peaked on this 
galaxy, may suggest that the group is in the late stages of dynamical 
relaxation. 

Emission is detected in the imaging data out to r = 8.9 arcmin 
(r ~ 230 kpc) at 5a above the X-ray background level evaluated 
from a surrounding annulus. A radial surface brightness profile is 
shown in Fig. f2] extracted from the centroid of emission in bins 
containing a signal-to-noise ratio of S/N 5. As expected from 
the irregular X-ray morphology, a standard /3-model is not a satis- 
factory description of this profile, with \l of 3.5 for 52 d.o.f. The 
data show significant deviations from the best-fitting model (with 
/3 = 0.39 ± 0.01 and r c = 21.5+5 J arcsec) at all radii. However, 
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the fit residuals do not exhibit any systematic radial variation, sug- 
gesting they are caused by local fluctuations in the IGM distribution 
rather than large-scale inhomogeneities. Hence, despite the fact that 
the /3-model fit is clearly statistically unacceptable, it remains use- 
ful for our purposes as a means of characterizing the global hot gas 
distribution. From inspection of Fig. [2] it is also not clear that a dif- 
ferent, or more complex, model would be able to provide a better 
description. 

The relatively broad point spread function (PSF) of XMM, not 
accounted for in the surface brightness fitting, could potentially af- 
fect the observed profile at small radii, and hence the derived core 
radius and central gas density. We do not expect this to be an im- 
portant effect, however, because even just the innermost radial bin 
in Fig.[2]extends to r = 12 arcsec, roughly twice the full-width at 
half maximum of the EPIC PSF. The fact that the best-fitting core 
radius is another factor of two larger also suggests that PSF blurring 
does not have a significant impact on the derived results. 

Using the double-subtraction approach for extracting a back- 
ground spectrum, a fit to the global 0.3-5 keV spectrum extracted 
inside r = 6 arcmin (r = 150 kpc) gives a temperature T — 
°- 83 ±o'.06 keV and abundance Z = 0.03 ± 0.01 Z , thus confirm- 
ing the low abundance derived fr om ROSAT data within the same 
region dOsmond & Ponmanl2 004). These values are consistent with 
those obtained using local background subtraction, but results are 
better constrained due to the superior statistics of the blank-sky 
background data. The derived flux and surface brightness profile 
imply a central hot gas density no = 4.4 ± 0.5 x 10~ 3 cm -3 . 

4.1.3 HCG30 

Despite this group representing the most Hi deficient system within 
our sample, the Chandra data do not reveal any clear evidence 
for diffuse IGM emission. No extended sources outside individual 
galaxies are detected by 'vtpdetect', thus corroborating the RASS- 
based result o f P96. With a galaxy velo city dispersion of only 
72 km s _1 , the Osm ond & Ponmanl d2004h scaling relations would 
suggest a very low IGM temperature ofT = 0.2±0.1 keV. In order 
to test for the presence of any such gas, we generated 0.2-0.4 keV 
images of the data on the S2 and S3 CCDs separately. These im- 
ages were exposure-corrected and smoothed but not background- 
subtracted, in an attempt to suppress any bias related to ACIS cal- 
ibration uncertainties at these low energies. The results reveal no 
clear spatial variations in the diffuse emission on either chip, sug- 
gesting emission at a level consistent with the local background. 

As a further test, we searched the unsmoothed data for a radial 
gradient in the exposure-corrected 0.3-2 keV emission level across 
the S2 and S3 CCD's, finding no significant variation with distance 
from the optical group centre (see Fig. [3] for a surface brightness 
profile extracted on the S3 CCD). This implies that any diffuse IGM 
emission would have to be near-uniformly distributed on scales of 
~ 300 kpc, an unlikely scenario for this low-a system, in which the 
angular extent of the region encompassing the four principal group 
members is only ~ 5 arcmin (~ 90 kpc). 

The RASS 0.5-0.9 keV count rate in a 0.5-1 deg annulus cen- 
tred on the optical group centre is 8a above the exposure-weighted 
mean of the appropriate Chandra blank-sky data, suggesting a con- 
siderable contribution from either background or (Galactic) fore- 
ground emission at this position. The high background rate relative 
to blank-sky data requires us to evaluate the background from the 
source data, thus restricting the source region under investigation 
on S3 to within r ~ 3 arcmin (~ 50 kpc) of the optical group cen- 
tre. Assuming that the background in the data can be safely evalu- 



ated from source-free regions on S3 outside this central region (an 
assumption supported by Fig. [3), the density of any hot IGM in the 
group can be constrained. Inside this region, and for T in the range 
0.1-0.3 keV and any subsolar metallicity, the 3a upper limit to the 
mean gas density is (n e ) < 1.2 x 10 -3 cm -3 . The possibility 
that T is very low in this group propagates into a relatively weak 
constraint on (n e ). 



4.1.4 HCG37 

Irregular diffuse X-ray emiss ion was detected in th is group with 
ROSAT out to r ~ 8 arcmin dMulchaev et al.ll2003h , well beyond 
the region covered by a single ACIS chip in our Chandra data. 
The latter clearly indicate that the group emission is sharply peaked 
on the early-type galaxy HCG 37a (= NGC 2783), the nucleus of 
which is also detected as a point-like source in the data. The as- 
sociation of the IGM X-ray peak with HCG 37a was not obvious 
from t he earlier ROSAT data, as the much broader ROSAT PSF re- 
quired iMulchaev et al. I d2003l) to exclude point-like emission out 
to r = 1.5 arcmin from the peak, thus effectively masking out 
HCG 37a in the data. Despite the overall irregularity of the group 
emission, the X-ray centroid (when masking out the HCG 37a nu- 
cleus) coincides to within 10 arcsec with the optical position of 
HCG 37a as listed in NED. 

Since group emission covers the S3 CCD, we cannot reliably 
use any method relying on source-free regions on S3 to evalu- 
ate the background in the Chandra data. The situation is aggra- 
vated by the fact that RASS data indicate a 3a soft background 
deficit at this position relative to the appropriate blank-sky data, 
so background subtraction by means of these is not straightfor- 
ward either. To circumvent these is sues, we adopted the method 
employed bv lVikhlinin et al.1 d2005b . First, a source minus blank- 
sky spectrum was extracted on the back-illuminated SI chip, to 
quantify the difference in the soft background between source- and 
blank-sky data. The spectrum was fitted with a T = 0.18 keV 
Z = Zq mekal plasma in the 0.4-1 keV range, with the nor- 
malization allowed to be negative. The best-fitting model was then 
added to the model fit of the blank-sky subtracted source emission 
on S3 inside r — 3.7 arcmin (r ~ 100 kpc) after scaling to the 
source region area. The resulting background level was also used 
for the surface brightness analysis. We note that, at 90 per cent 
confidence, the best-fitting T and Z resulting from this approach, 
T = 0.86±°;Jg keV and Z = O.W±q° 8 4 Z q , are just consistent 
with the P96 values inside r = 150 kpc (T = 0.67 ± 0.11 keV, 
Z — 0.17 ± 0.15 Zq), lending some credibility to this approach. 

The surface brightness profile shown in Fig. [2] confirms the 
presence of emission across the full S3 CCD in the Chandra data. 
The profile was centred on the X-ray peak and extracted in bins 
containing at least 30 net counts. Despite the irregularity of the 
emission on large scales (Fig. QJ, a /3-model provides a good fit 
to the profile across the full radial range plotted in Fig. [2] yielding 
/3 = 0.42 ± 0.02 and r c = 3.5±{ l arcsec, with xl = 0.80 for 
19 d.o.f. The spatial and spectral results imply a gas density in the 
group core of 0.038 ± 0.006 cm" 3 . 



4.1.5 HCG40 

The Chandra observation of this group was split into two sepa- 
rate pointings, so a merged event file was produced for the imaging 
analysis. Although undetected in a 3.6-ks ROSAT pointing (P96), 
Fig.Q]suggests the presence of diffuse emission in this group. The 
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centroid of this emission, with the optical extent of the individual 
group members masked out, is located ~ 0.5 arcmin to the NW 
of HCG40c and so is not clearly associated with any individual 
galaxy. To test that the emission seen in Fig.Q]is truly extended and 
not simply due to the smoothing of point sources, a surface bright- 
ness profile of the unsmoothed, exposure-corrected emission was 
extracted from the centroid in bins of at least 30 net counts, with 
individual galaxies masked out. The result is shown in Fig. [2] with 
the background evaluated from off-source regions on the S3 CCD. 
Emission is detected above this background out to r — 2.3 arcmin 
(r = 65 kpc), suggesting group-scale extended emission, although 
the detection is only significant at > 3<j for the innermost 30 kpc. A 
/3-model provides an acceptable fit to this profile, with xt = 1-01 
for 7 d.o.f., yielding f3 = 0.60io. 4 ° and r c = 46.9±f|;? arcsec 
(23lg 7 kpc), in accordance with expectations for a typical X-ray 
bright group. 

For the spectral analysis of this emission, spectra and response 
products were extracted separately for each of the two observations. 
The spectra were then jointly fitted within the central r = 1 arcmin, 
within which the signal allows useful constraints to be obtained, 
using a surrounding 2.5-3.5 arcmin annulus for background esti- 
mation. With only ~ 140 net counts, the IGM abundance remains 
unconstrained. Fixing Z at 0.3 solar yields T — 0.59loli2> and T 
remains consistent with this for any subsolar Z. For these param- 
eters, the observed flux translates into a central electron density of 
l.l±0.2x 10~ 3 cm -3 and implies a diffuse 0.3-2 keV luminosity 
inside r = 2.3 arcmin of 3.1 ± 0.5 x 10 40 erg s _1 . 

The extent of the emission, coupled with the fact that it is 
not clearly centred on any group member, suggests that the emis- 
sion is not due to, for example, hot gas associated with an el- 
liptical but rather reflects the presence of a hot IGM. This in- 
terpretation would place HCG40 among the relatively rare ex- 
amples of spiral-dominated groups showing intergalactic hot gas; 
within Hickson's (1982) catalogue, only HCG 16, 57, and the well- 
studied HCG92 (Stephan's Quintet) share similar features (e.g., 
Dos Santos & Mamonll999l : lFukazawa et alj2002l ; lTrinchieri et all 
2003b . Based on the B-band luminosities of the group mem- 
bers, and on the Lx- Lb relations for ellipticals and n orma l star- 
forming spirals f rom | p' Sullivan. Ponman & Collins! d2003l) and 
iRead & PonmarJ J200lb~ respectively, one would expect a total 
galactic diffuse Lx ~ 9 x 10 40 erg s _1 in the group, a factor 
of three larger than that found here for the intragroup emission. Al- 
though care has been taken in masking out emission from the group 
members, the low S/N and the compactness of the galaxy config- 
uration implies that we cannot exclude a residual contribution to 
the diffuse emission from individual galaxies. A conservative ap- 
proach would be to regard the association of the observed diffuse 
emission with an intragroup medium in HCG 40 as tentative rather 
than conclusive. 

4.1.6 HCG 44 

While Fig. Q] does not indicate the presence of any IGM emission 
in this system, this could simply be an artefact of the proximity of 
the group (D ~ 23 Mpc) in combination with the limited Chandra 
angular coverage, which furthermore renders quantitative analysis 
of the background level in the data non-trivial. Using local back- 
ground subtraction could potentially produce unreliable results, as 
IGM emission might cover the entire ACIS array. The situation is 
further complicated by an enhanced particle level in the cleaned 
data compared to the blank-sky files, with the 10-12 keV count 
rate on the S3 CCD being 35 per cent higher than the correspond- 



ing blank-sky value. In addition, soft Galactic 0.5-0.9 keV emis- 
sion at a level of 2a above the blank-sky data is also present at 
this position, so it is not obvious that blank-sky data would be an 
appropriate choice for background estimates. 

Fortunately, the presence of an overlapping XMM pointing 
allows an independent test for the presence of diffuse emission 
in the group. The Chandra and XMM surface brightness profiles 
shown in Fig. [3] suggest no detectable IGM emission close to the 
optical group centre. Note that the extraction of these profiles ex- 
cluded different position angles due to the optical group centre be- 
ing close to the southern (northern) edge of the S3 (EPIC) CCDs, 
so the profiles extend in largely opposite directions on the sky. The 
XMM profile shows no systematic variation out to r — 15 arcmin 
(r ss 100 kpc), remaining largely consistent with the background 
level evaluated outside this region in the source data. Consequently, 
the background level for the Chandra profile was estimated from 
the northern corners of the S3 CCD, the result suggesting no ex- 
cess diffuse emission extending northwards either. 

Furthermore, no extended sources that can be unambiguously 
associated with group emission were detected by 'vtpdetect'. In ad- 
dition to the two group galaxies HCG 44a and b, a third extended 
X-ray source is seen on the S3 chip, clearly visible in Fig.Q]roughly 
four arcmin north of the spiral HCG 44a, and also seen in both the 
XMM data and in pointed ROSAT observations. The X-ray peak of 
this source coincides with a 2MASS source with vtik = 14.31, 
but there is no optical counterpart or redshift information available 
in NED. If this source were at the group distance, the resulting 
LT-band luminosity of 2.1 x 10 8 Lq,k would place it at the ex- 
treme faint end of the dwarf galaxy luminosity function, with the 
ROSAT flux implying a ratio Lx / Lk ~ 0.013, two orders of mag- 
nitude above typical values seen even for dwarf starburst galaxies 
( IRasmussen. Stevens & Pon man 2004). Spectral fit results provide 
further support for the idea that this source is unlikely to be associ- 
ated with HCG 44. A thermal plasma model fixed at the group red- 
shift returns an unacceptable fit (xt = 1.61), whereas a significant 
fit improvement results when leaving z as a free parameter, yield- 
ing X l = 1-37 for 6 d.o.f, with T = 1.98t° 6 3 l and z = °- 24 ^o'.07 
for an assumed abundance of Z = 0.3 Zq. 

While the Chandra data are useful in terms of investigating 
evidence for hot gas being stripped from individual galaxies in this 
very nearby system, it is not clear that these data enable significant 
improvements on the hot IGM constraints over the existing 4.7- 
ks ROSAT pointing (with its much larger field of view enabling a 
more reliable background subtraction in this case). Even the XMM 
data can only probe emission from a quarter of the volume inside 
r — 100 kpc from the optical group centre, due to the latter being 
close to the edge of the XMM field of view. Using the adopted cr-T 
relation, which would suggest T = 0.4 ± 0.1 keV, the ROSAT 
constraint of P96 on the X-ray luminosity inside r = 150 kpc 
(Lx < 3.6 x 10 40 erg s _1 for our adopted distance) translates into 
(n e ) < 1.0 x 10~ 4 cm -3 for any subsolar metallicity. The cor- 
responding XMM constraint of (n e ) < 1.4 x 10~ 4 cm~ 3 applies 
within a volume less than 10 per cent of that probed by ROSAT, so 
we have adopted the stronger ROSAT limit in Table[2] 

4.1.7 HCG 97 

This constitutes the most X-ray luminous system within the sam- 
ple. A two-dimensional ana lysis of a 21-ks XMM observation of 
this group was performed by Mahdavi et al. (2005), along with op- 
tical spectroscopy identifying 37 members. Their XMM data show 
a plume stretching to the southeast, beyond the region covered by 
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Fig.Q] Mah davi et al.l {2005 ) speculate that this plume could repre- 
sent gas stripped from one of the central galaxies, but the XMM data 
alone cannot establish this, and the Chandra data on the S2 CCD 
cannot improve on the situation beyond confirming the presence 
and overall morphology of this feature. One of the spectroscopi- 
cally identified member galaxies is located within the plume but is 
not itself detected in either the XMM or Chandra data. On the S3 
CCD, the emission appears fairly regular in Fig. [T] albeit with the 
central X-ray contours somewhat elongated towards the south-east. 
If masking out this elongated feature, the centroid of the IGM emis- 
sion coincides to within 10 arcsec with the position of the optically 
brightest group galaxy, HCG 97a, as listed in NED. 

As established already by ROSAT observations (e,g., 
iMulchaev etafl [2003). diffuse emission in this group extends well 
beyond the region covered by the S3 CCD, so blank-sky data were 
used to evaluate the background for the Chandra surface bright- 
ness analysis (the RASS 0.5-0.9 keV background count rate at this 
position is in good agreement with the corresponding exposure- 
weighted mean value of the blank-sky data). We note, however, that 
the Chandra observation was somewhat affected by background 
flares, reducing the useful exposure time from 57.9 to 36.3-ks. As 
for HCG 44, the background remains high after cleaning, with the 
10-12 keV count rate on S3 again being 35 per cent above the 
blank-sky value. While bearing this issue in mind, results indicate 
that emission is detected above 5<r significance everywhere on the 
S3 chip. A fit to the exposure-corrected 0.3-2 keV surface bright- 
ness profile, extracted from the X-ray peak and shown in Fig. [2] 
yields /3 = 0.414 ± 0.004 and r c = 6.6 ± 0.6 arcsec, with 
xt = 1.45. However, while Figure [2] and the exc ellent agreement 
of the se results with the best-fitting parameters of IMulchaev et al.l 
d2003l) (who find (3 = 0.41 ± 0.01 and r c < 0.1 arcmin) sug- 
gest that our background estimate is not seriously in error, we will 
nevertheless base our normalization of the density profile on the 
ROSAT results of IMulchaev et all J2003l) . given the concern about 
the elevated particle background in the Chandra data. Combin- 
ing their spectral results and X-ray luminosity (0.3-2 keV Lx = 
1.20^0 24 x 10 42 er 8 s_1 f° r our adopted distance) with our sur- 
face brightness fit then implies a temperature T = 0.97^']^ keV 
and central density of 0.048 ± 0.002 cm -3 , which are the values 
listed in Table[2] 

4.1.8 HCG 100 

This is a group in which the Hi is clearly being stripped from the 
galaxies at present, with much of it pulled into a 100 kpc long 
tidal tail extending to the southwest from the optical group centre 
(Borthakur et al., in prep.). In addition, the VLA data reveal a strik- 
ing Hi trail extending to the east away from the group core, pro- 
truding from one of the galaxi es in the field, Mrk935. This galaxy 
is not included in the original HicksqjJ dl982l) catalogue, but is a 
group member on the basis of its projected distance from the optical 
group centre (6.7 arcmin ~ 130 kpc) and small radial velocity dif- 
ference of Av ~ 250 km s" 1 relative to the group mean, as listed 
in NED. The associated Hi feature may therefore indicate ongoing 
stripping as the galaxy falls into the group. HCG 100 thus consti- 
tutes an excellent laboratory for the study of the processes whereby 
Hi is removed from individual galaxies and heated. Unlike the case 
for the other Chandra observations presented here, this group was 
observed using the ACIS-I array, to allow the observed field to fully 
encompass all of the interesting Hi features mentioned above. 

The Chandra observation was split into two, so the imaging 
analysis proceeded as for HCG 40. The combined imaging data, 



shown in Fig. [TJ and the resulting surface brightness profile in 
Fig.[3] do not reveal any clear indications of diffuse emission above 
the background level as evaluated outside r — 8 arcmin from the 
corners of the ACIS-I array. Although Fig. [3] indicates a mild net 
excess in the two innermost bins, the signal within this region is 
significant at less than 1.5a. We also note that no extended sources 
are detected outside individual galaxies with 'vtpdetect', and that 
the group also remained undetected in RASS data (P96). With the 
a~T relation suggesting T = 0.3 ± 0.1 keV, the absence of an 
IGM detection inside r ~ 100 kpc (r m 5 arcmin) from the optical 
group centre implies a 3u upper limit on the mean IGM density of 
(n e ) < 3.3 x 10" 4 cm -3 . 

4.2 Individual galaxies 

The results presented so far demonstrate the presence of a de- 
tectable hot IGM within half of our sample only. However, even in 
the absence of such gas, there could still be an intragroup medium 
present with temperature or density below our detection limits (in- 
cluding any Hi already stripped from individual galaxies, as evi- 
denced by the GBT detection of intergalactic Hi in many of our 
groups). 

To explore the possibility that group galaxies could be inter- 
acting with such a medium, and to search for signs of galaxies be- 
ing stripped of any hot gas, we present in Fig. [4] a collage of all 
group members which were clearly identified as X-ray extended 
sources by our source detection algorithms. These images were 
adaptively smoothed following the procedure outlined in Section[3] 
Exceptions are HCG 30a and Mrk 935, for which a simple Gaussian 
smoothing (with a — 10 arcsec) was employed due to the very low 
S/N. For the groups observed by XMM (HCG 7 and 15), the combi- 
nation of group distance and the broader EPIC PSF does not enable 
a clear distinction between point-like and diffuse emission, so none 
of the relevant galaxies has been included in this figure. HCG 44b, 
lying on a chip gap in the Chandra data and not covered by the 
overlapping XMM pointing, has also been excluded. 

The figure does not reveal any clear evidence for galaxies cur- 
rently being stripped of any hot gas. In particular, there are no in- 
dications of X-ray tails or bow-shock features indicating interac- 
tions with a surrounding medium. Such ta ils have been observed 
extending from galaxies in c lusters (e.g.. IWang. Owen & Ledlowl 
|2004 ISun & Vikhlininll2005T) but seem to be very rare in groups, 
with pe rhaps NGC 6872 and NGC 2276 the most pr ominent ex- 
amples dMachacek et alj2005l ; lRasmussen et alj20 06). There is an 
indication of an asymmetric structure in HCG 97d even in the un- 
smoofhed data, with a hint of a tail pointing south, but the sig- 
nal is too weak to exclude contamination by a faint point source. 
Mrk 935, with the remarkable Hi tail extending to the east, also 
presents evidence of some X-ray asymmetry in this direction, but 
the S/N is again too low to enable firm conclusions. 

In addition to diffuse galactic emission, X-ray point sources 
in individual galaxies are detected in all groups, with some of these 
sources clearly associated with the galaxy nuclei. The relevance 
of this is linked to the possibility that galaxies suffering strong 
(tidal) interactions could be showing enhanced nuclear activity, for 
example associated with a nuclear starburst or strong AGN accre- 
tion fueled by a tidally induced gas inflow. For reference, Table [3] 
lists detected X-ray sources whose position coincides with the op- 
tical centre of individual group members, with luminosities of any 
point-like component derived as described in Section [3] Note that 
in some cases, such as HCG 97e and most of the XMM sources, we 
cannot clearly distinguish between nuclear and galaxy-wide diffuse 
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Figure 4. Smoothed 0.3-2 keV images of individual group galaxies show- 
ing diffuse X-ray emission. A horizontal bar marks a scale of 1 arcmin in 
each case. 



emission, and the classification of these is followed by a '?' in the 
Table[3] For HCG 7, however, the tentative identification of nuclear 
X-ray activity in three of the four princ ipal galaxies agrees per fectly 
with the Spitzer far-infrared results of Gall agher et al.l ( I2008I) . sug- 
gesting our identification is reasonably robust. 

Among the principal members in Hickson's (1982) catalogue, 
we detect 22 candidate nuclear X-ray sources in 40 galaxies, with 
roughly two-thirds of these falling in the groups with a detectable 
IGM. The corresponding nuclear source fractions of 64 ± 17 and 
44 ± 16 per cent in the groups with and without detectable hot 
gas, respectively, are statistically indistinguishable at the la level. 
If instead splitting the sample according to Hi deficiency, the cor- 
responding fractions are 42 ± 15 (high Am) and 67 ± 18 per cent, 
a difference which is only just significant at la. The median nu- 
clear X-ray luminosities for the two subsamples are also very sim- 



Table 3. Overview of X-ray sources centred on individual group galaxies 
as identified by our detection algorithms. Galaxy morphologies were taken 
from NED. Column 4 lists the unabsorbed 0.3-2 keV luminosity of any 
nuclear component. 
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ilar, 4.4 X 10 39 (high A ffl ) and 4.9 x 10 39 erg s" 1 , suggesting 
that the above conclusions are not strongly biased by a system- 
atic difference in limiting X-ray flux between high- and low-Am 
groups. We also note that results from optical spectroscopy indi- 
cate that ~ 40 per cent of the principal members in HCGs in gen- 
eral show evidence for AGN activity, with a total of ~ 70 per cent 
showing emission line s from either AGN or star formation activity 
(Martinez et al. I l2007h . Our derived fractions are generally brack- 
eted by these values, suggesting that our results give a reasonably 
reliable picture of the frequency of nuclear activity within our sam- 
ple. 

With the limited statistics available, there is thus no strong 
evidence from the X-ray data alone for enhanced nuclear activity 
within a certain kind of groups in our sample. Specifically, if inter- 
preting the X-ray bright or highly Hi deficient systems as dynami- 
cally more evolved, we find no clear indication that the frequency or 
strength of nuclear X-ray activity depends on the dynamical status 
of the group. However, we note that this result applies to a small 
sample and to the principal members only; a complete census of 
galaxy membership from optical spectroscopy would be required 
to extend this analysis to optically fainter group members and place 
this conclusion on a more robust basis. The tentative lack of a clear 
enhancement in nuclear X-ray activity among the most Hi defi- 
cient groups within our sample may tie in with the observation that 
star formation activity is not globally enhanced in HCG galaxies 
compared to isolated ones ( Verdes-Montene gro et al.ll9 98). as dis- 
cussed in more detail in Section|6] 



12 J. Rasmussen et al. 



4.3 Hi deficiency and IGM properties 

The observed diversity in the diffuse X-ray properties of these 
highly Hi deficient groups immediately suggests that galaxy-IGM 
interactions are not the dominant mechanism for driving cold gas 
out of the galaxies within our sample and establishing the observed 
Hi deficiencies. Of course, this conclusion neglects the fact that we 
are not uniformly sensitive to the presence of hot gas in the differ- 
ent groups. A quantitative comparison of observed Hi deficiencies 
and derived IGM properties is therefore presented in Figure|5] The 
left panel shows Ahi and hot IGM mass as listed in Table[2] with 
both quantities derived within the same region (inside rm in the 
Table). Even when considering the X-ray detected systems alone, 
the obvious lack of a positive correlation between the two quanti- 
ties immediately suggests that the amount of hot gas in the group 
core is not a pivotal factor for Hi removal. We note that an identical 
conclusion is reached if replacing Ahi with 'missing' Hi mass in 
the plot. The strong 3cr upper limit (M < 2.6 x 10 9 M Q ) on the 
IGM mass in the highly Hi deficient HCG 7 (with the largest 'miss- 
ing' Hi mass in the sample of ~ 1.8 x 10 10 Mq) only reinforces 
this conclusion. Also note, as pointed out in Section |4~Tl that more 
realistic assumptions about the IGM distribution in the X-ray unde- 
tected groups could reduce their upper limits to A/igm by perhaps 
an order of magnitude, but that this has no bearing on the above 
conclusions. 

Hence, neither is it surprising that Ahi does not show a 
clear dependence on the characteristic IGM ram pressure plotted 
in Fig. [5p and evaluated as the product of ay and the volume- 
weighted mean IGM density Migm/V within the volume V = 
(4/3)7rr HI covered by the radio data, with all quantities taken from 
Table [2] Note that the very compact HCG 40 - in which we cannot 
completely rule out a residual galactic contribution to the derived 
IGM mass - stands out among the X-ray detected groups, with 
a characteristic ram pressure 1-2 orders of magnitude below that 
seen in the X-ray bright systems. For the X-ray undetected groups, 
similar comments apply as for Fig. [5^. 

Finally, in Fig. [5J: we investigate the dependence of Ahi 
on IGM temperature. Thermal evaporation of galactic Hi through 
heat conduction from the IGM is expected to proceed at a rate 
M oc T 5 / 2 if unsuppressed by, for example, magnetic fields. Given 
this strong temperature dependence, the lack of a positive Ahi-T 
correlation suggests that heat conduction is not an important effect 
within our sample. The location of the exceptionally Hi deficient 
HCG 30 in Fig. [5J; would seem to pose a particular challenge for 
this mechanism. 

Overall, Fig. [5] therefore seems to confirm the notion that 
Hi deficiency is not tightly linked to the presence or nature of an 
IGM in these groups. There are some caveats to this interpreta- 
tion though. For example, it is worth emphasizing that it is the four 
groups with the highest velocity dispersion that are X-ray detected. 
If the remaining groups contain warm (T < 10 6 K) rather than 
hot gas, and thus fall well below the a-T relation for X-ray bright 
systems, our constraints on the hot gas density could seriously un- 
derestimate the true IGM density in these systems. Unfortunately, 
this possibility cannot be directly tested with the present data. How- 
ever, the fa ct that our X-ray de t ected systems scatter fairly tightly 
around the lOsmond & Ponm an ( 2004) relation, as shown in the in- 
set in Fig. [5};, may support our use of this relation for predicting 
T also for the undetected systems. Furthermore, results for other 
X-ray detected groups suggest, if anything, that the poorest sys- 
tems tend to have T on the high side for their velocity dispersion 



dOsmond & Ponmani r2004) although the situation could, of course, 
be different for groups that remain X-ray undetected. 

Another concern is that cry an d hence P ram in Fig. |5p may 
not be robustly determined, being based on just a handful of bright 
galaxies in most cases. Large corrections to oy would be needed, 
however, in order to affect our overall conclusions. A further point 
is that any intergalactic Hi already stripped from individual galax- 
ies could potentially contribute to the IGM mass and ram pressure. 
However, for the X-ray detected systems, the total Hi mass inside 
the GBT beam is of order 5 per cent of the corresponding hot gas 
mass (Borthakur et al., in prep.), suggesting that any cold gas can 
be neglected for the present purposes. 

Despite the appearance of Fig. [5] it is premature to exclude 
the possibility that galaxy-IGM interactions could play a role for 
Hi removal in some of our groups. For example, ram pressure strip- 
ping is expected to occur when the IGM ram pressure exceeds the 
gravitational restoring pressure of the galaxy. The efficiency of this 
process therefore depends not only on the properties of the IGM, 
but also on those of the individ ual group galaxies. The process of 
viscous stripping ( TNulserJll982h shares similar features, and could 
be operating even when ram pressure itself is insufficient to remove 
any Hi. Finally, ram pressure may also indirectly affect the Hi in 
the disc. For example, many disc galaxy formation models pre- 
dict that massive spirals are surrounded by hot gaseous haloes from 
which gas may cool out to provide fuel for ongoing star forma- 
tion in the disc (e.g., lToftetalj|2 002). The removal of this coronal 
gas by external forces could contribute to Hi deficiency, if the lim- 
ited supply of Hi in the disc is consumed by star formation with- 
out being replenished from the hot halo (strangulation; see, e.g., 
iKawata & Mulch aev 2008). We next seek to quantify the impor- 
tance of these various mechanisms. 



5 THE IMPACT OF THE IGM: MODELLING 
GALAXY-IGM INTERACTIONS 

In an attempt to constrain the role of galaxy-IGM interactions for 
typical disc galaxies in the individual groups, we constructed sim- 
ple analytical models of galaxies moving through the hot intragroup 
gas in the gravitational potential of each group. As explained in 
Section |4~T1 the derived Hi deficiency for each group should be 
viewed as an average for the group members, since it is often non- 
trivial to evaluate observed Hi masses for the individual members. 
For our modelling purposes, we have therefore adopted a single, 
fiducial galaxy model, with overall properties broadly matched to 
the fairly well-constrained mean properties of the late-type group 
members in our sample. 

As described in detail below, the group potential and galaxy 
orbits are less well determined for each group. Consequently, we 
evolve the adopted galaxy model according to four different as- 
sumptions for each group, corresponding to two choices for the 
gravitational potential, and two for the galaxy orbits within the cho- 
sen potential. The variation among the resulting mass losses from 
the galaxy can serve as a means to gauge the uncertainties associ- 
ated with these assumptions. We will consider two different strip- 
ping processes for the cold gas in the disc, viz- classical ram pres- 
sure stripping and turbulent viscous stripping. For the gas in any 
hot halo, we only consider ram pressure stripping for simplicity. 
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Figure 5. Hi deficiency and (a) hot gas mass inside the region used for determining Ahi, (b) characteristic ram pressure, and (c) hot gas temperature for the 
va rious groups. Empty circles represent groups with no detectable hot gas. Inset in (c) shows velocity dispersion vs. Tjqm for the X-ray detected groups, with 
the lOsmond & Ponrn an 12004) relation, i.e. equation {3), overplotted as a dashed line. 



5.1 Group model and galaxy orbits 

Our goal is to evaluate the efficiency of ram pressure stripping 
and related processes, without resorting to detailed numerical mod- 
elling, which is beyond the scope of this work. For this pur- 
pose, we would ideally adopt a single value of the ram pressure 
for each grou p. The simplest app roach is to assume the classi- 
cal analytical iGunn & Got 3 dl972h stripping criterion along with 
a constant ram pressure equal to its peak value. Hydrodynami- 
cal simulations invo lving a constant ram pressure have shown the 
IGunn & Gottl dl972h criterion to be remarkably accurate in terms 
of predicting the disc strippin g radius and the mass of gas lost 
dAbadi. Moore & Bowed U999T) . and it remains a reasonable ap- 
proxi mation even when allowing for orbital variation s in ram pres- 
sure dJachvm et alj200llRo"ediger & Briiggenll2007l) . 

However, re cent hydrodynamica l simulations of disc galaxies 
moving in radial iJachym e t alj|2007l) and two-dimensional orbits 
jRoediger & Brii ggen 2007) within a non-uniform gas distribution 
have revealed an important exception to this rule. If the ram pres- 
sure changes faster than the characteristic stripping time-scale, as 
will be th e case for galaxies m oving through a highly concentrated 
IGM, the lGunn & Gottl d 19721) criterion tends to overestimates the 
stripping efficiency. Since stripping is not instantaneous, an ISM 
element may not always be accelerated to galactic escape velocity 
before the peak of the ram pressure is over. In such cases the gas 
wil l eventually re-accre te, a possibility not taken into account by 
the lGunn & Gottl d 1972b criterion. In the present study, these con- 
siderations could potentially be relevant for several of our groups, 
given the fairly small core radii resulting from the surface bright- 
ness fits. Instead of using the peak value of the ram pressure in 
our modelling, it therefore seems more sensible to adopt an orbit- 
averaged mean ram pressure by evaluating the time spent by the 
galaxy at a given velocity and IGM density. 

To this end, we first derived total mass profiles M to t(r) for 
each X-ray detected group, assuming a spherically symmetric gas 
distribution in hydrostatic equilibrium, 

»r / x -kT(r)r /dlnn(r) dlnT(r)\ 
Mt„t < r) = v ' - y - ' + - v ' . (4) 
Gfj,rn p \ amr amr J 

Since in general we have neither the statistics nor the spatial cov- 
erage to constrain T(r) to large radii, we made two assumptions 
about the temperature profile which are likely to bracket the actual 



temperature distribution in these somewhat disturbed groups. We 
either simply assumed T(r) equal to a constant mean value (T), or 
T(r) — — 0.5(T) log(r/r 50 o) + 0.67, app ropriate for reasonably 
undis turbed groups outside any cool core dRasmussen & Ponmanl 
1 20071) . In the latter case, rsoo, the radius enclosing a mean den- 
sity of 500 times the critical value p c , was evaluated iteratively un- 
til convergence was reached. For (T) we used the measured value 
listed in Table[2] 

The resulting mass profiles were cha racterized analytically by 
fitting them with standard 'NFW models dNavarro. Frenk & White! 
1 1997 ), in which the dark matter distribution is described by 

p( r ) = Pc 7-7 w., ^ , s 2 , 0) 

[r/r s )(l + r/r s y 

where 8 C is a dimensionless parameter related to total group mass 
M200, and r s is a scale radius which reflects t he more commonly 
used concentration parameter c = T200 / r s (see lNavarro et aljl997l 
for details). With A/200, and hence S c , fixed from the measured 
(extrapolated) mass profile, equation l[5) was fitted to this profile 
to derive values of c under the two assumptions about T(r) de- 
scribed above. The results, summarized in Table [4] imply typical 
values of M200 » 1 — 2 X 10 13 M Q with c ~ 5 - 10. For these 
group masses, derived concentration parameters are thus in good 
agree ment with expectatio ns from cosmological TV-body simula- 
tions dBullock et al.ll200"lh . Note that the assumption of a declin- 
ing temperature profile typically reduces the derived group mass 
by ~ 30 per cent while increasing the halo concentration by a fac- 
tor of ~ 2. 

For the galaxy orbital configuration within the derived gravi- 
tational NFW potential, we assume two different orbits, both radial 
and hence going through the group core. The galaxy is assumed 
to be experiencing a face-on IGM encounter in either case. These 
maximizing assumptions allow us to estimate how important ram 
pressure stripping can ideally be in our groups. The two orbits dif- 
fer only in terms of the assumed initial position and velocity of the 
galaxy, with the galaxy initially at rest at a small clustercentric ra- 
dius in the first scenario, and falling towards the group core from a 
large radius and with a high initial velocity in the second. Specifi- 
cally, the following two scenarios are considered: 

(i) For each group, we determine f, the observed (projected) 
mean clustercentric distance of the principal galaxies in each group. 
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Table 4. Results of NFW fits to the derived group mass profiles under the 
two assumptions for 7'(r) described in the text. 
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The model galaxy is assumed to be initially at rest at a larger clus- 
tercentric distance rrj, from which it falls freely towards the group 
centre, ro is chosen such that when the galaxy reaches r = f, it has 
attained a velocity corresponding to the observed group velocity 
dispersion. Typical values are ro ~ 100 kpc and f ~ 40 kpc. 
(ii) The galaxy enters the group halo at ro = r2oo with a radial 
velocity v r corresponding to the halo circular velocity at this radius, 
v r = (GM /r2oo) 1 ^ 2 ■ Typical values are ro ~ 500 kpc and v r ~ 
400kms _1 . 

In both cases we follow the galaxy until it turns around, having 
completed one passage through the group core. The first scenario is 
chosen in an effort to match the observed average position (modulo 
projection effects) and galaxy velocity in each group at present. In 
practice, it represents galaxy motion fairly close to the group core, 
with a mildly varying ram pressure, and so is somewhat reminis- 
cent of the 'classical' ram pressure scenario involving a constant, 
high ram pressure. The assumption underlying this orbit is extreme, 
however, in the sense that the true clustercentric distances will gen- 
erally be larger than the observed (projected) ones, which implies 
that galaxies will generally spend a larger fraction of their time at 
large distance than implied by this assumption. Therefore we also 
consider scenario (ii) as a kind of opposite extreme. In this, galaxies 
experience a considerably higher peak ram pressure, but this occurs 
only relatively briefly. In the following, these two scenarios will be 
referred to as orbit (i) and (ii), respectively. 

Although these two orbits do not necessarily encompass the 
extreme orbital solutions for the group members, they do represent 
two rather different cases, thus offering a handle on the uncertainty 
in the predicted mass loss related to orbital assumptions. Further- 
more, while completely radial orbits may not be very common, we 
note that Hi deficient spirals in clusters tend to have more eccentric 
orbits than non-deficient ones dSolanes et alJuOOll) . and that cos- 
mological infall along filaments would proceed in fairly eccentric 
orbits, thus lending some support to our simplifying assumption. 
For a detailed discussion of the impact of orb ital parameters on the 
stripping efficiency, we refer to lHesteJ feOOrj) . 



5.2 Galaxy model 

For the galaxy model, needed to estimate the gravitational restor- 
ing force and hot halo thermal pressure of a 'typical' late-type 
group galax y within our samp l e, we follow the general approach 
described in iRasmussen et al.l 12006) which is repeated here for 
completeness. The model consists of a spherical dark matter (DM) 
halo with density profile 



Ph(r) = 



M h n exp(-r 2 /n 2 ) 



27r 3/2 nr 2 (i + r»/r;|) ' 
a hot gaseous halo of the same form, a spherical bulge with 

M b 1 



Pb(r) = 



27rr^ r(l + r/rb) 3 ' 



and exponential stellar and gaseous discs, each of the form 

M d 



Pd(R,z) 



exp(— _R/i? d )sech {z/z d ) 



(6) 



(7) 



(8) 



Here Mb, Mh, and Md are the total masses of each component, rb 
and rh are the scalelengths of the bulge and halo, respectively, r t is 
the DM halo 'truncation' radius, R d is the cylindrical scalelength 
of the disc components and z d the corresponding thickness, and 
finally 



r, = {1 - ^ 1/2 gexp(g 2 )[l - erf(g)]}- 1 , 



(9) 



where q — r\/r t and erf is the error function. With this model, the 
restoring gravitational acceleration §j(-R, z) in the direction z per- 
pendicular to the disc can be eval uated analytically f or each model 
component using the equations of lAbadi et al.l a 19991) . to whom we 
refer for more details. 

In order to constrain model parameters, stellar masses of the 
HCG members were evaluated from their B- and JC-band mag- 
nitudes as listed in N ED, following the prescription adopted by 
iMannucci et al.l d2005h . For this purpo se, only the prin cipal mem- 
bers in our eight groups, as listed by Hickson ( 1982), with mor- 
phological types later than SO were included, yielding a mean stel- 
lar mass of 4.1 x 10 1() Mq. For a subset of these galaxies (8 
out of 27), maximum disc rotational velocities, indicative of total 
galaxy masses, are also available in the Hyperleda database, with 
a mean value of 135 km s" 1 . We note that, in terms of mean stel- 
lar mass, this subset is representative of the full sample, showing 
(M») = 4.4 x 10 10 M Q . We therefore assume a stellar mass of 
4 x 10 10 Mq for the galaxy model, distributed such as to yield a 
bulge-to-disc mass ratio of 1/4, appropriate for an Sb/c spiral. Using 
the relation of lHavnes & Giovanellil ( fl984h . we further assume an 
Hi mass of 6.7 x 10 9 Mq to ensure that our model galaxy initially 
has a 'normal' Hi mass for its sample-averaged blue luminosity of 
L B = 1.6 x 10 10 L . 

The scalelengths of the stellar and gaseous disc components 
in the model were chosen to ensure that at least 85 per cent of the 
stellar and Hi mass resides within the average optical disc radius of 
rr> — 10 kpc, as derived from the size of the D25 ellipse for each 
spiral member in our groups. For simplicity, the gas distribution in 
the hot gaseous halo is assumed to follow that of the underlying 
dark matter, but with a smaller value of r t corresponding to twice 
the 'size' r_o of the stellar disc, and with a total mass corresponding 
to the Hi mass in the disc. The hot halo is assumed to be isothermal 
at a temperature T ~ 0.06 keV, the virial temperature correspond- 
ing to the maximum allowed disc rotational velocity in the model, 
taken to be ~ 150 km s _1 . Once the baryonic model components 
have been specified, the parameters of the DM halo are effectively 
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Table 5. Adopted parameters of the galaxy model. L is the characteristic 
scalelength for each component, i.e. r h and rt for the haloes, rj, for the 
bulge, and R d and z d for the disc components. 



Component A^total L 

(10 10 Mq) (kpc) 



DM halo 


28 


5, 100 


Hot gas halo 


0.7 


5, 20 


Stellar bulge 


0.8 


0.5 


Stellar disc 


3.2 


3, 0.25 


Cold gas disc 


0.67 


3, 0.25 



set by this maximum allowed rotation velocity, and by the require- 
ment that the model baryon fraction match the universal value of 
~ 15 per cent. Table [5] summarises the adopted model galaxy pa- 
rameters. Note that the cold gas component in the model refers to 
the distribution of Hi only, as we do not consider any molecular 
gas here. In summary, this model roughly reproduces the average 
stellar mass, bulge-to-disc ratio, disc size, and maximum disc rota- 
tional velocity seen for the spirals in our groups, with a total baryon 
fraction consistent with the universal value, and an initial Hi mass 
as expected for a non-stripped isolated galaxy with these properties. 



5.3 Computing mass losses 

Combining equation l[5} with the measured density distribution of 
intragroup gas, time-averaged values of IGM density (p) and the 
square of the orbital velocity (v 2 al ) experienced by the galaxy in 
its orbit can be evaluated. For the averaging time-scale, we only 
consider the segment of the orbit for which the ram pressure ex- 
erts a significant influence on the ISM, taken to be from the time at 
which the instantaneous ram pressure, if sustained, would remove 
at least 1 per cent of the cold gas. This is to avoid the artificial sup- 
pression of (p) and {v 2 ) that would otherwise result from including 
the time spent by the galaxy at large radius where ram pressure is 
completely negligible. Note that the corresponding characteristic 
ram pressure is independent of group mass and orbital initial con- 
ditions, as it depends only on the assumed galaxy model. 

Assuming a face-on IGM encounter, the condition for ram- 
pressure stripping is then evaluated as 

*{^ + ^ + -k + -k) < v>m> (10) 

where ^(R,z) is the restoring gravitational acceleration in the 
direction z perpendicular to the disc, originating from the stellar 
bulge (subscript 'b'), dark matter halo ('h'), gaseous disc ('g') and 
stellar disc ('*'), respectively, and E g is the s urface density of cold 
gas. Equation dlOt is similar to the classical iGunn & Gottl d 19721) 
stripping criterion, but takes into account the mass distribution in 
the galaxy rather than simply assuming a homogeneous disc thin 
enough to be described solely by its surface density. Its solution 
also provides us with the 'stripping region', the region in the (R, 
z)-plane from which gas is permanently lost by the galaxy. 

Apart from conventional ram pr essure s tripping, transport 
processes such as viscous stripping dNulsenl 1 1983) . caused by 
Kelvin-Helmholtz instabilities arising at the ISM-IGM interface, 
could also play a role even when th e ram pressure itself is in- 
sufficient to remove g alactic gas dOuilis. M oore & Bowed l200d : 
iRasmussen et alj[2006l) . Turbulent viscous stripping of a gas disc 



of radius ro is expected to operate at Reynolds numbers Re > 30, 
where 

Re = Mr D /\, (11) 

M is the Mach number of the IGM flow past the galaxy, and A oc 
T 2 n _1 is the ion mean free path in the IG M. The expected mass- 
loss rate due to this process dNulserJI 19821) . 

M vs ~ 0.57rr-opw ga i, (12) 

scales only linearly with galaxy velocity and so could be important 
in a wider range of environments than ram pressure itself. We in- 
clude this process in the stripping calculations by evaluating equa- 
tions dl It and d!2l l at each point in the orbit and adding up the 
total mass loss. Dealing wi th disc galaxie s rather than spheroids, 
we have used only half the Nulser] Tl982l) mass loss rate in equa- 
tion d!2t because of the correspondingly smaller galaxy surface 
area for a given r — ro- For ro itself, we use the smaller of 
the sample-averaged value of ro = 10 kpc and the stripping ra- 
dius predicted by equation dlOl l. Note that M ^ 1 in equation dl lb 
even if the galaxy is moving supersonically, as the post-shock IGM 
flow past the galaxy will always be subsonic. Also note that the 
stripping efficiency of both ram pressure and viscous stripping 
should be largely unaf fected by the presence of a shock front (see 
IRasmussen et al j|2006l) . 

Having dealt with the stripping of cold gas from the disc, we 
now turn to the removal of any gas situated in a hot galactic halo. 
It seems plausible that, at the very least, the cooling and gradual 
inflow of any such gas will be disrupted once the external ram 
pr essure P ram exceeds its thermal pressure P t h- The simulations 
of Mori & Burkert] j200Ct) confirm that this condition provides a 
reasonable estimate of the mass of a galaxy that will have its hot 
halo completely stripped by ram pressure. For the purpose of also 
assessing the importance of strangulation for the model galaxy, we 
therefore use the derived values of (p) and {v 2 } to evaluate the 
'strangulation region' where P ram > P t h, and the corresponding 
mass of affected coronal gas. The aim here is only to develop a 
rough picture of the potential importance of IGM interactions for 
strangulation, so the effects of viscous stripping are not considered 
for the hot halo. 

We readily acknowledge that our modelling approach is infe- 
rior to detailed numerical models (e.g., Hester 2006), and hydrody- 
namical simulations in particular. It is our hope that one can nev- 
ertheless have some confidence in the results, given that we have 
employed a reasonably detailed galaxy model and are making some 
allowance for orbital variations in ram pressure. The advantage of 
the adopted method is that it is easily tailored to the specific condi- 
tions in individual groups at little computational cost. 

5.4 Results 

We re-iterate that we are considering four scenarios for each group, 
consisting of two separate assumptions about the galaxy orbit, and 
two about the IGM temperature (and hence total mass) distribution 
in the groups, as specified in Section 15.11 The results of the Hi 
stripping calculations for each of these four cases are summarized 
in Table[6] which lists the derived r.m.s. orbital velocity w rm s, the 
orbit-averaged IGM number density, the mass of Hi lost due to ram 
pressure and viscous stripping, and the total fraction of Hi stripped 
from the initial model reservoir of 6.7 x 10 9 Mq. Recall that v TjnB 
is computed only from the point in the orbit where ram pressure 
becomes significant and so depends not only on orbital parameters 
but also on the IGM distribution. 
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Table 6. Orbit-averaged galaxy velocities and IGM densities, along with 
predicted Hi mass loss AAf due to ram pressure ('rp') and viscous stripping 
('vs') after one passage through the group core for the assumed orbits and 
group mass profiles. The final column lists the total fraction / of Hi lost. 



Orbit ( i), T isothermal 



Group 


^rms 


(n) 


AM rp 


AM va 


/ 




(kms" 1 ) 


(cm 3 ) 


(10 9 M Q ) 


(10 9 M s ) 




HCG15 


354 


6.7 x 10" 4 


0.9 


1.0 


0.29 


HCG37 


342 


1.5 x 10~ 3 


1.2 


1.4 


0.39 


HCG40 


199 


7.8 x 10" 4 


0.6 


0.2 


0.12 


HCG97 


362 


2.6 x lO- 3 


1.6 


2.3 


0.58 


Orbit ( i), T declining 


Group 


^rms 


<n> 


AM rp 


AM va 


/ 




(kms" 1 ) 


(cm" 3 ) 


(10 9 M Q ) 


(10 9 M s ) 




HCG15 


404 


7.2 X 10" 4 


1.0 


0.9 


0.29 


HCG 37 


402 


1.9 x 10" 3 


1.4 


1.2 


0.39 


HCG40 


266 


8.2 x 10" 4 


0.8 


0.2 


0.14 


HCG 97 


435 


2.8 x 10~ 3 


1.8 


2.0 


0.57 


Orbit ( ii), T isothermal 


Group 


^rms 


(n) 


AM rp 


AM va 


/ 




(kms" 1 ) 


(cm 3 ) 


(10 9 M Q ) 


(10 9 M ) 




HCG 15 


498 


1.2 X 10" 4 


0.6 


1.2 


0.27 


HCG 37 


541 


1.5 x 10~ 4 


0.7 


1.9 


0.39 


HCG 40 


630 


4.3 x 10~ 5 


0.5 


0.3 


0.12 


HCG 97 


536 


2.8 x 10" 4 


0.9 


4.3 


0.77 


Orbit ( ii), T declining 


Group 


"rms 


(n) 


AM rp 


AM VS 


/ 




(kms" 1 ) 


(cm 3 ) 


(10 9 M Q ) 


(10 9 M ) 




HCG 15 


448 


1.2 X 10" 4 


0.5 


1.2 


0.26 


HCG 37 


491 


1.5 x 10" 4 


0.6 


1.9 


0.38 


HCG 40 


554 


4.6 x 10~ 5 


0.4 


0.3 


0.11 


HCG 97 


477 


2.8 x 10" 4 


0.8 


4.3 


0.76 



As can be seen from the Table, the amount of gas lost through 
either stripping process is generally an appreciable fraction of the 
initial Hi mass. Viscous stripping in particular can remove a sub- 
stantial fraction of the cold disc gas for all model assumptions. 
HCG 40 is the exception, with the relatively tenuous IGM in this 
group removing at most 10-15 per cent of the Hi. The steeper IGM 
density profile in this group also implies that a larger fraction of the 
total IGM mass is encountered at high velocity, leading to a higher 
AA:frp/AM vs ratio than for the other groups, and a higher value 
of Mrms in orbit (ii), because the ram pressure becomes significant 
closer to the core in this system. 

The Table further shows that the main variation in the com- 
puted mass loss for a given group and stripping mechanism derives 
from the choice of orbit rather than the assumed group mass pro- 
file. Orbit (i) is generally more efficient at removing gas through 
ram pressure stripping than orbit (ii), because although the peak 
ram pressure is considerably higher in the latter case, the galaxy 
spends a comparatively shorter time in regions corresponding to 
high values of Pram - Conversely, viscous stripping is more efficient 
in orbit (ii), as this mechanism acts even at relatively low v ga ,i, al- 
lowing the associated mass loss to build up significantly over the 
much longer crossing time-scale relevant for this orbit. Note, how- 
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Figure 6. Isobars of gravitational restoring pressure (dashed) and hot halo 
thermal pressure (dotted) for our fiducial model galaxy in the different 
groups. Dashed lines outline the galactic regions outside which Hi in the 
disc can be stripped by ram pressure, corresponding to equality in equa- 
tion Mot . Dotted lines show the corresponding regions for the hot halo gas. 
The order of the contours is the same in both cases. 

ever, that for a given group, the outcome in terms of total Hi mass 
loss, AJVftot = AM r p + AM VS , is largely insensitive to the various 
orbit and mass profile assumptions, despite clear variations in (n) 
and «rms among the four scenarios. The only exception to this is 
the significantly more massive HCG 97, for which the deeper grav- 
itational potential and higher IGM mass implies that orbit (ii) is 
relatively more efficient at removing galactic gas than for the other 
groups. 

Figure [6] outlines the stripping region for a model galaxy in 
each of the X-ray detected groups, based on solving equation d 10b . 
The figure shows the one of our four scenarios in which the effect of 
ram pressure is generally most pronounced, i.e., orbit (i) with T(r) 
declining. In the outer disc, the gravitational restoring pressure, and 
hence the stripping region for the cold gas, is seen to be nearly in- 
dependent of vertical disc height for interesting values of \z\, and 
the gravitational restoring force at a given R peaks well above the 
disc. Compared to the 'size' of the stellar disc (ro = 10 kpc), it 
is clear that the gas disc becomes mildly truncated by ram pressure 
in HCG 97 and 37 but remains largely unaffected in the other two 
groups. Note that this truncation reduces the viscous stripping effi- 
ciency by reducing the surface area of the gas disc exposed to the 
IGM. Viscous stripping is therefore slightly more efficient for low 
values of P ra m, adding to the explanation of the higher values of 
AM VS for orbit (ii) in Table[6] 

Predicted Hi deficiencies corresponding to the mass losses in 
Table [6] are compared to the observed values in Fig. [7] Errors on 
predicted deficiencies correspond to the full range in predicted Hi 
mass loss for each group under the various orbit and mass pro- 
file assumptions. Shown are the expectations from ram pressure 
stripping alone, as well as the full Hi mass loss from combining 
equations J10l > and J12b . The dashed line in the figure represents 
equality between modelled and observed Hi deficiencies; under our 
model assumptions, anything below this line cannot be explained 
by galaxy-IGM interactions alone. The plot clearly suggests that 
ram pressure stripping on its own is not sufficient to cause the ob- 
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Figure 7. Observed Hi deficiencies compared to the model predictions of 
Table [6] for the stripping of Hi by ram pressure alone (empty diamonds), 
and by ram pressure plus viscous stripping (shaded). Dashed line represents 
equality between observed and predicted Am- 



Table 7. Fraction of hot halo gas lost due to ram pressure in the various 
orbital scenarios. 



Group 


Orbit (i) 


Orbit (i) 


Orbit (ii) 


Orbit (ii) 




T(r) const. 


T(r) decl. 


T(r) const 


T(r) decl. 


HCG15 


0.53 


0.64 


0.31 


0.27 


HCG37 


0.74 


0.85 


0.39 


0.35 


HCG40 


0.31 


0.46 


0.21 


0.18 


HCG97 


0.84 


0.94 


0.52 


0.48 



served Hi deficiencies, even in this X-ray bright subsample of our 
groups. When including viscous stripping, however, if as efficient 
in removing Hi as assumed here, galaxy-IGM interactions can cer- 
tainly help explain observed values of Ahi, but they can potentially 
fully account for the Hi loss only in HCG 37 and HCG 97, i.e. in 
just two out of our eight groups. 

Regarding the issue of strangulation, Fig.[6]also illustrates the 
derived stripping region for hot halo gas according to the adopted 
fram > Pth criterion. The figure suggests that ram pressure alone 
could remove a sizable fraction of the halo gas in the model. Ta- 
ble |7J lists the fractions of stripped halo gas in the different scenar- 
ios, showing that these can be substantial for orbit (i) in particular, 
peaking at ~ 95 per cent in HCG 97. Viscous stripping, not in- 
cluded here, could potentially contribute beyond these estimates. 
This suggests that galaxy-IGM interactions in X-ray bright groups 
could play an important role in removing the gas supply that may 
ultimately fuel star formatio n in spirals, in qualitative a greement 
with the simulation results o f lKawata & Mulchaevl |2003) . For our 
specific model setup, there is a large dispersion in the fraction of 
gas affected, however, and it is not clear that the effect would be 
important in groups such as HCG 40. 



5.5 Model limitations and caveats 

It is important to stress that the calculations presented here are only 
intended to provide a rough picture of the impact of galaxy-IGM 
interactions in our sample, and we do not claim that these results are 
anything but indicative. Caveats include the fact that the adopted 
galaxy model parameters provide a plausible, but not necessarily 
unique, representation of the spirals in our sample. It is also a 
possibility, albeit one we cannot easily evaluate, that some group 
members may be individually less Hi deficient than the 'group- 
averaged' value of Ahi, and so could have Hi deficiencies con- 
sistent with removal by ICM interactions alone, even if our model 
results suggest otherwise for the group members as a whole. 

As regards the model calculations, an important limitation is 
the fact that the model is completely static. In practice, the gas 
distributions both in the groups and their galaxies will be evolv- 
ing, which could be particularly relevant for orbit (ii) where orbital 
time-scales are several Gyr. Also, the estimate of viscous stripping 
mass loss does not take into account that this process itself will re- 
duce the size of the gas disc and hence the mass loss rate according 
to equation J 1 2b . Furthermore, this stripping process may saturate 
jNulsen|[l982f7 which would reduce the mass loss, the presence of 
magnetic fields could also suppress hydrodynamical instabilities, 
and the presence of a hot halo could perhaps to some extent shield 
the cold gas disc from such instabilities. The latter could be of par- 
ticular relevance for orbit (ii), in which AM VS is high but where a 
comparatively smaller fraction of halo gas is lost due to ram pres- 
sure. These considerations suggest that the estimated contribution 
by viscous stripping to the Hi mass loss should be regarded as an 
upper limit. 

With their compact galaxy configurations and low velocity 
dispersions, our groups could also represent environments in which 
galactic dark matter haloes are subject to significant tidal trunca- 
tion. Hence, another issue is to what extent our results are affected 
by the adopted assumptions on dark matter in the galaxy model. 
For the adopted disc and bulge parameters, our freedom to modify 
the assumed DM distribution is mainly constrained by the require- 
ment that the maximum disc rotational velocity u max of the model 
should not exceed the allowed ~ 150 km -1 . A full exploration of 
model parame t er spa ce is beyond the scope of this work, and we 
refer to Hester) (2006) for a more thorough investigation of these is- 
sues. Nevertheless, to provide a rough picture of the impact of DM 
for our results, we repeated the stripping calculations of Section l53l 
with three different modifications, subject to the above ii m „ crite- 
rion: 

(i) Assuming no DM [or larger DM halo scalelengths rt and rn 
in equation ([6)]. This must be regarded as an extreme assumption, 
which should facilitate stripping. 

(ii) Increasing the DM mass by a factor of two (for which the v max 
criterion then requires at least one of the DM scalelengths rt and 
Th to increase by a similar factor). This should suppress stripping. 

(iii) Assuming a higher DM concentration (i.e. lower ru or ri) by 
a factor of three (suppresses stripping), which in turn requires a 
lower total DM mass by a factor of five (facilitates stripping). These 
choices comply with the typical tidal truncation of DM ha los in- 
ferred for galaxies in massive clusters I lLimousin e t al. 2007). 

We find that, in all three cases, the amount of Hi lost by the model 
galaxy is at most modified by 10-15 per cent, regardless of the 
group and orbit considered. This perhaps somewhat surprising re- 
sult has its origin in the fact that the vast majority of cold gas 
in the model is located at low galactocentric distances, where the 
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restoring gravitational acceleration - which determines the strip- 
ping region in accordance with equation l l 10b - is dominated by 
the baryonic disc components and not the dark matter halo. The 
variations in the stripped gas mass for the hot halo can be larger, 
up to 45 per cent for HCG 40, but they can generally easily be ac- 
commodated by the variations associated with the different orbital 
assumptions. We therefore conclude that our results are reasonably 
robust to changes in the assumed amount and distribution of dark 
matter in the galaxy model. 



6 DISCUSSION 

The groups in our sample are all Hi deficient even when account- 
ing for any intergalactic Hi not clearly associated with individual 
galaxies. Removal of Hi from the group members is by itself insuf- 
ficient to explain this situation, as the removed gas must also be pre- 
vented from staying neutral. Thus, unless Hi is somehow destroyed 
in situ within the galaxies, it must go through a two-stage process 
whereby it is first removed from the galaxies and then ionized by 
a possibly unrelated mechanism. Here we discuss these different 
possibilities in light of our X-ray and modelling results. 

6.1 In situ destruction of Hi 

The Hi could potentially have been destroyed within the galaxies 
themselves, either by evaporation though thermal conduction from 
the IGM, or through consumption by star formation, provided there 
is no continuous replenishment of cold material. A third possible 
mechanism involves heating and possibly ejection from the galax- 
ies by starburst winds. While the first possibility, direct heating by 
the IGM, seems implausible in light of the absence of a positive cor- 
relation between Am and Tigm (cf. our discussion of Fig.[5j;), the 
two other scenarios deserve some further attention. Hi consumption 
by star formation could help establish observed deficiencies, but a 
prerequisite is that the consumed gas is not simply replenished at 
a similar rate, for example through the cooling out of hot, coro- 
nal material on to the disc. Our model results for the X-ray bright 
groups (Table[7J suggest that such a strangulation scenario could at 
least be greatly facilitated by the removal of coronal gas due to ram 
pressure, but a more detailed exploration of model parameter space 
would be required to assess the general validity of this conclusion. 

In this context, it is instructive to compare the results 
from our simple an a lytical model to the simulation results of 
iKawata & Mulchaevl (2008). These authors investigate the effi- 
ciency of ram pressure stripping and strangulation for a galaxy 
in a small galaxy group on the basis of a cosmological smoothed 
particle hydrodynamics (SPH) simulation. They do not specifically 
consider a 'compact' group, but this has the advantage that their 
target galaxy is not subject to noticeable gas loss from tidal interac- 
tions (D. Kawata, priv. comm.), so in this sense a comparison to our 
model calculations is justified. Their target galaxy has initial prop- 
erties broadly similar to those of our galaxy model in terms of stel- 
lar, gaseous, and total mass. It is followed for one passage through 
the group as it enters a group of virial mass M ~ 8 x 10 12 Mq 
from an initial position roughly corresponding to twice the virial 
radius. The galaxy orbit thus shows some similarities to our orbit 
(ii), but, with a pericentre at r ~ 200 kpc, do es not take the galaxy 
throug h the very core of the IGM distribution. IKawata & Mulchaevl 
(2008) find that the resulting ram pressure induces strangulation to 
a significant degree. Star formation and thus Hi consumption be- 
comes mildly enhanced during infall, and the hot halo gas is almost 



completely removed during the group passage. The combination of 
these processes has effectively consumed the Hi and quenched star 
formation by the time the galaxy re-emerges at the virial radius. 

Our model conclusion that ram pressure stripping may 
strongly affect any coronal gas in the X-ray bright groups is thus in 
encouraging agreement with these simulation results, and it is even 
possible that our model underestimates the importance of ram pres- 
sure in this context. This is made more relevant still if ram pressure 
- or other galaxy interactions with the group environment - indi- 
rectly accelerate strangulation by enhancing the disc star formation 
rate and hence the consumption of Hi. This possibility draws ob- 
servational support from the strong star formatio n activity seen in 
the ra m-pressure affected group spiral NGC 2276 (Rasmussen et al. 
2006), as well as from recent hydro-simulations dKronberger et al.l 
2008). However, Verd es-Montenegro et alj ( fl998l) noted on the ba- 
sis of star formation rates derived from IRAS far-infrared luminosi- 
ties that there is no indication of enhanced star formation among 
the HCG galaxies compared to the level seen in isolated galaxies. 
In fact, current star formation rates in these groups generally seem 
too low to explain the missing Hi by consumption through star for- 
mation. For the spirals in our sample, a mean star formation rate of 
M* « 1.4 M0 yr _1 can be derived for the /iM5-detected galaxies. 
However, many of our galaxies remain undetected in one or more 
IRAS bands and so have only upper limits to M« (if including these 
upper limits in the mean, the result is M* < 0.9 M© yr _1 ). At 
these rates, the time-scale for star formation to exhaust an initial Hi 
supply of 6.7 x 10 9 Mq to current levels is at least 5 Gyr, even if 
neglecting the return of unprocessed material to the ISM and any 
cosmological accretion of gas by the galaxy. 

Thus, while o u r mo del calculations and the results of 
IKawata &~M ulchaev (2008) suggest that hot halo gas can be re- 
moved fairly efficiently in several of our groups, potentially inhibit- 
ing or at least suppressing any replenishment of disc Hi from a hot 
halo, strangulation does not seem to have played an important role 
in establishing current Hi levels within our sample. At the observed 
star formation rates, the remaining Hi can continue to fuel star for- 
mation for many Gyr without drastically reducing the Hi supply. 
Of course, star formation rates could have been much higher in 
the past, which could also indirectly have contributed to exhausting 
the gas supply through the ejection of gas from the disc by star- 
burst winds. However, unless M» in these galaxies has generally 
declined to current levels only fairly recently or has been affected 
by the group environment over cosmological time-scales, the impli- 
cation seems to be that destruction of neutral hydrogen within the 
galaxies themselves cannot explain the Hi deficiencies. The pres- 
ence of intergalactic Hi in some of our groups also shows that this 
scenario cannot provide an exhaustive explanation. We are there- 
fore compelled to also consider the alternative, externally-driven 
removal and destruction of Hi. 

6.2 Hi removal by external forces 

Focussing first on the removal of Hi, both galaxy-IGM and galaxy- 
galaxy interactions could be envisaged to play a role. Our results 
demonstrate that although earlier studies indicated a link between 
significant Hi defi ciency in Hickson groups and the presence of 
hot intragroup gas JVerdes-M ontenegr o et alj200ll) . there is clearly 
no one-to-one correspondence between the two. The lack of a 
detectable intragroup medium inferred here for half of the eight 
most Hi deficient compact groups specifically appears to rule out 
IGM dynamical interactions as generally dominant for Hi removal 
within these environments. Ram pressure and viscous stripping 
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could nevertheless still have played a role for Hi removal in our X- 
ray detected systems. Our modelling results indicate, however, that 
the efficiency of these processes is generally insufficient to fully ac- 
count for the Hi missing from the individual group members, espe- 
cially if our simplistic treatment of viscous stripping overestimates 
the predicted Hi loss, as hypothesized in Section [53] 

Overall, the simulations of Ka wata & Mulchaevl J2008h seem 
to support these conclusions. In these, ram pressure affects the 
amount of cold gas in the disc to an even lesser degree than in 
our model, although this could perhaps be attributed to the fact 
that our groups are at least twice as massive, and perhaps also to 
the maximizing orbital assumptions adopted in our model. Viscous 
stripping d oes not appear to be impor tant in the simulations, but as 
noted by Kawata & Mulchaev (2008), such processes are not nec- 
essarily well treated by SPH schemes, so a direct comparison to our 
results may not be meaningful. 

The limited impact of galaxy-IGM interactions inferred for 
our groups commands an alternative explanation for the Hi re- 
moval. Tidal stripping constitutes an obvious candidate in these 
dense, low-cr environments, particularly in light of the fact that 
the four X-ray undetected groups within our sample exhibit the 
lowest velocity dispersions and so could be expected to represent 
environments where tidal interactions should be most important. 
However, even if for now ignoring the problem of subsequently 
heating the removed Hi, it is not immediately clear that tidal in- 
teractions would necessarily result in increased Hi deficiency, as 
they would also affect the stellar component in the galaxies. If stars 
and gas are removed in roughly equal proportion, and if the re- 
moved stars become part of any unde tected intracluster light (cf. 
iGonzalez. Zaritskv & Zabludoal2007h while the Hi remains de- 
tectable, this could potentially even result in an Hi excess. Recall- 
ing that Hi deficiency is defined here on the basis of the observed 
B-band galaxy luminosity, increasing Ahi through tidal stripping 
therefore requires either the preferential removal of cold gas com- 
pared to stars (and its subsequent heating), or that Lb is simultane- 
ously boosted relative to the Hi mass for the non-stripped compo- 
nents. 

The latter possibility gains support from the observation that 
specific star formation rates generally tend to be higher in galaxies 
with close neighbours JLi et alj |2008). with enhanced nuclear star 
formation plausibly arising as a consequence of tidally induced gas 
inflow. It is perhaps curious then, as mentioned above, that there is 
no evidence for enhanced star formation in Hickson groups com- 
pared to the level in isolated galaxies. This runs contrary to expec- 
tations for interacting galaxies, and would seem to argue against 
Lb being significantly boosted in the HCG members relative to 
their remaining stellar or Hi mass. The other possibility mentioned 
above, that Hi is more easily tidally stripped than the stellar com- 
ponent, is perhaps more promising. In many of our groups, such as 
HCG 100 (cf. Section B. 1.8t . significant amounts of intergalactic Hi 
is detected with the GBT, whereas the stellar components are not 
noticeably affected. A possible explanation is that a relatively larger 
fraction of Hi compared to stars initially resided at large galacto- 
centric radii, where tidal removal would be most efficient. The fact 
that the Hi disc in typical spirals and late-type dwarfs is often more 
extended than its optical counterpart, with the radial distribution of 
Hi declining less steeply with radius than that of the B-band light 
jBroeil s & Rhee 1 997 ; ls waters et al]|20 02), seems to support such 
an explanation. 

Tidal stripping clearly is taking place in some of our groups, 
notably in the X-ray undetected groups HCG 100 and HCG 44. In 
the latter, both optical (Fig. [TJ and Hi data (Borthakur et al., in 



prep.) show strong evidence for the SBc galaxy HCG44d being 
tidally stripped. This may indicate that tidal stripping is the primary 
mechanism by which Hi is removed from the galaxies in these two 
groups. Although it is tempting to extend this conclusion to all the 
X-ray undetected groups in our sample, and perhaps even beyond, 
the inconspicuous star formation rates in Hickson groups in gen- 
eral, and the tentative lack of enhanced nuclear X-ray activity in 
the highly Hi deficient systems in particular (Section[4j2}, may not 
argue in favour of such an explanation. 



6.3 Destruction of removed Hi 

As emphasized in the beginning of this Section, Hi deficiencies 
can only be explained if the hydrogen, once removed from the 
group members, is also transformed from its neutral phase. Irre- 
spective of the processes accomplishing its removal, a mechanism 
must therefore also be invoked for ionizing the Hi during or follow- 
ing its transfer to intergalactic space. In our X-ray bright groups, a 
candidate process is readily available, since any removed Hi is ex- 
pected to evaporate due to heating by the ambient IGM at a rate 
M oc r 5 / 2 , almost i ndependently of realistic values of the IGM 
density JSpitzenl 19621 ). A detailed comparison of the X-ray and Hi 
properties of these groups should help test this explanation and will 
be presented elsewhere (Verdes-Montenegro et al., in prep.). In the 
X-ray undetected groups, where any IGM is expected to be rela- 
tively cool (cf. Fig. |5j;) if at all present, the picture is less clear- 
cut. Even if these groups do contain an IGM, the lower predicted 
IGM temperatures in this subsample imply average IGM heating 
time-scales an order of magnitude above those in the X-ray bright 
systems, casting doubt on whether this mechanism would be suffi- 
cient. 

If gas is predominantly removed by tidal interactions in the 
X-ray undetected groups, it is instead conceivable that some of the 
Hi has been heated by tidal shocks, although the presence of inter- 
galactic Hi in some of the groups would imply heating time-scales 
well in excess of those associated with the Hi removal itself. An- 
other possibility is that the column density of any removed Hi be- 
comes too low for the gas to be self-shielding against ionization by 
cosmic UV radiation. If so, much of the undetected hydrogen in 
these groups could potentially be in the form of a tenuous, photo- 
ionized intergalactic plasma. A quantitative investigation of these 
possibilities is the subject of future work, but at present a clear pic- 
ture of the fate of the removed Hi in the X-ray undetected groups 
remains elusive. 



7 CONCLUSIONS 

Based on a sample of eight Hickson compact groups selected 
for their high Hi deficiencies, we have used Chandra and XMM- 
Newton data to assess the properties of any hot intragroup medium 
(IGM) and constrain the role of galaxy-IGM interactions in remov- 
ing Hi from the galaxies in these groups. The X-ray analysis re- 
veals a detectable IGM in four of the eight groups. We have tenta- 
tively identified the detected diffuse emission in HCG 40, a spiral- 
dominated group, as associated with an intragroup medium, but 
the combination of a low signal-to-noise ratio and an exceptionally 
compact galaxy configuration precludes a highly robust conclusion 
for this particular system. The remaining three groups are all fairly 
X-ray luminous, showing substantial amounts of intergalactic hot 
gas with a somewhat disturbed morphology in all cases. 
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The remarkable X-ray diversity seen across the sample imme- 
diately suggests that the presence of a significant IGM is not a dom- 
inant factor in establishing observed Hi deficiencies, despite ear- 
lier re sults indicating such a connection ( Verdes-Montenegr o et al.l 
l200lh . It is particularly notable that some of the most Hi deficient 
groups show no detectable hot IGM, including HCG 30 which only 
contains a few per cent of the expected Hi mass for its galaxy con- 
tent. A comparison of Hi deficiency with either hot IGM mass or 
characteristic 'mean' ram pressure confirms the lack of a clear cor- 
relation even for the X-ray bright systems, although statistics are 
naturally too limited to enable firm conclusions on this basis alone. 
The Hi deficiency does not seem to depend on IGM temperature ei- 
ther, suggesting that heat conduction from the IGM does not play an 
important role in destroying galactic Hi (although once removed, 
this gas is likely to evaporate on fairly short time-scales in the X- 
ray bright groups). 

From fitting analytical models to the derived mass profiles of 
the X-ray detected groups, we have constructed plausible models 
of the gravitational potential and associated radial galaxy orbits for 
each of these groups. Combined with the inferred IGM distributions 
and a numerical model of a late-type galaxy with properties broadly 
matching those of our observed spirals, this has enabled estimates 
of the importance of ram pressure stripping and viscous stripping in 
removing Hi from the late-type galaxies in each group. The results 
indicate that, even under maximizing assumptions about the galaxy 
orbit, ram pressure stripping will remove only small amounts of 
cold gas from the group members, peaking at 10-25 per cent in the 
X-ray bright HCG 97. We find that viscous stripping is generally 
more efficient, with the combination of the two processes capable 
of removing more than half of the cold ISM in HCG 97 and po- 
tentially fully accounting for the missing Hi mass in both HCG 37 
and HCG 97. However, the efficiency of viscous stripping is likely 
overestimated with our simple analytical approach, and yet these 
processes are insufficient in terms of explaining the Hi deficiency 
of the X-ray detected HCG 15 and HCG 40. 

The model results also indicate that ram pressure can effi- 
ciently remove a large fraction of any hot galactic halo gas that 
may otherwise act as a supply of fresh material for star forma- 
tion in the disc. However, even if the gas supply to the disc can 
be completely cut off, gas consumption at the typical star forma- 
tion rates in the groups would proceed far too slowly to explain 
the observed shortfall of Hi by itself. Much higher star formation 
rates in the recent past are required for this process to have had any 
significant impact. This may suggest that the observed Hi deficien- 
cies are not caused by in situ destruction of Hi within the galaxies 
themselves. It remains a possibility that even modest star forma- 
tion activity could have heated some of the Hi and lifted it above 
the disc midplane where it would be more susceptible to removal 
by ram pressure, simila r to the situation proposed for NGC2276 
(Rasm ussen et al.ll2006l) . The absence of observational signatures 
of this process, e.g., in the form of a hot gas tail extending from 
any of the group members, may suggest that such a mechanism is 
not generally very important within our sample though. 

By the process of elimination, it seems plausible that tidal in- 
teractions have played a key role for Hi removal in the groups, par- 
ticularly in those systems containing no detectable IGM. In order 
to explain the Hi deficiencies, this scenario would likely require 
preferential removal of Hi over stars, as perhaps facilitated by the 
more extended distribution of cold gas relative to stars in typical 
late-type galaxies. While it is perhaps not surprising that tidal in- 
teractions are affecting the gas content of galaxies in these compact 
groups, the tidal stripping explanation still faces some outstanding 



issues. Among these is the expectation that such interactions would 
generate enhanced star formation or nuclear activity, but there is 
no indication that the X-ray faint or highly Hi deficient systems in 
our sample show evidence for increased such activity (although the 
frequency of nuclear acti vity in galaxies in com pact groups in gen- 
eral may be rather high; Martinez et al. 2007). If interpreting the 
X-ray bright or highly Hi deficient systems within our sample as 
dynamically more evolved, we thus find no clear evidence that the 
frequency or strength of nuclear X-ray activity in the group mem- 
bers depends on the dynamical status of the group. It also remains 
unclear whether tidal interactions themselves can destroy the Hi 
during or following its removal from galactic discs and so fully ex- 
plain the Hi deficiency in any of our groups. 

In closing, our results suggest that galaxy-IGM interactions 
can have played a role for the removal and destruction of Hi in 
some of our groups, but a complete understanding of the origin 
of the observed Hi deficiencies and the processes causing it is 
still lacking. Strangulation or thermal evaporation do not emerge 
as important contenders, and typical indirect signatures of tidal 
interactions, such as enhanced star formation or nuclear X-ray 
activity, are not more pronounced within the more Hi deficient 
half of our sample. The latter s eems in line with previous results 
dVerdes-Montenegro et al.l 19 98) which indicate that star formation 
rates in Hickson compact groups are not globally enhanced relative 
to the field. We note here that this result could potentially be mis- 
leading, however, perhaps masking an evolutionary trend in which 
galaxies initially experience enhanced star formation which is then 
followed by an environment-driven suppression. A detailed corre- 
lation of Hi and X-ray morphology in the groups, coupled with a 
broad comparison of individual galaxy properties such as specific 
star formation rates, may therefore shed further light on the fate 
of the missing Hi in these compact systems. This is the subject of 
future work. 
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